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ABSTRACT 


A  finite  element  code  is  developed  for  analysis  and 
design  of  three  dimensional  truss  and  frame  structures. 
Structures  are  designed  for  minimum  weight  subject  to 
constraints  on:  member  stresses,  Euler  buckling,  shell 
buckling,  joint  displacements  and  system  natural  fre¬ 
quencies  .  Structures  are  optimized  with  respect  to  member 
size  and  structure  configuration. 

The  finite  element  code  may  be  used  for  analysis  only, 
or  may  be  coupled  to  an  optimizer  of  the  user's  choice. 

The  displacement  method  is  used  for  static  analysis ,  and 
structure  natural  frequencies  are  calculated  via  the  sub¬ 
space  iteration  method. 

Flexibility  is  provided  for  expansion  to  other  than 
tubular  frame  elements,  and  provisions  are  made  for  the 
future  growth  to  panel  and  other  types  of  structural 
elements . 

Documentation  is  provided  to  facilitate  use  of  the 
code.  A  User's  manual  is  presented  with  examples  and  re¬ 
sults.  An  explanation  of  how  this  code  may  be  coupled  to 
an  optimizer  is  also  provided. 
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The  task  of  the  Engineer  is  to  provide  the  best  solution 
to  the  customer’s  problem.  Usually  the  ’’best"  solution  is 
the  one  that  does  the  job  with  an  adequate  margin  of  safety, 
is  aesthetically  pleasing,  and  is  economically  feasible. 

The  solution  may  be  reached  through  various  means ,  but  effi¬ 
cient  use  of  design  tools,  as  well  as  efficient  use  of 
materials  is  important  since  both  add  to  the  overall  cost  of 
the  product.  Design  optimization  is  one  method  that  can  be 
used  to  maximize  the  efficiency  of  a  structure  by  minimizing 
its  weight  and,  presumably,  cost. 

Optimization  of  structures  has  had  continuing  changes 
since  its  development  in  the  early  1960 's  with  an  active 
area  of  research  being  elastic  structures.  The  main  goal  is 
to  design  structural  systems  that  efficiently  perform  speci¬ 
fied  purposes.  Since  most  physical  structures  can  be 
modeled  by  the  finite  element  method,  a  computer  program 
can  be  written  to  perform  the  necessary  calculations  to 
solve  the  problem. 

The  purpose  of  this  research  was  to  develop  a  finite 
element  code  that  could  be  used  to  analyze  a  combined  truss/ 
frame  structure  and  could  be  easily  coupled  to  an  optimizer; 
thus  providing  a  useful  tool  for  designing  such  structures 
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as,  for  example,  ships'  masts.  This  code  expands  the  pre¬ 
vious  work  by  Fitzgerald  [Ref.  1]  on  truss  structures  to 
the  more  general  six  degree-of-freedom  truss/frame  case. 

The  design  problem  considered  in  this  study  is  the  opti¬ 
mization  of  three-dimensional  statically  indeterminant  com¬ 
bined  truss-frame  structures  under  varying  load  conditions 
and  subject  to  a  variety  of  constraints.  The  objective  is 
to  minimize  the  weight  of  the  structure  where  the  design 
variables  are  member  sizes  and  joint  coordinates.  Constraints 
include:  maximum  normal  stress;  maximum  joint  displacements; 

minimum  structure  natural  frequencies;  Euler  buckling;  and 
in  the  case  of  the  tubular  frame  elements ,  local  or  shell 
buckling. 

In  the  present  code,  all  gradient  information  is  calcu¬ 
lated  by  the  finite  difference  method.  Modification  of  the 
code  to  permit  calculation  of  gradients  analytically  has 
been  identified  as  a  necessary  future  extension. 

This  document  describes  the  use  and  capabilities  of  the 
finite  element  computer  code  to  be  coupled  to  an  optimizer. 

The  user's  manual  presented  in  Chapter  V  contains  a  simple 
design  example  in  which  the  program  is  coupled  to  the  CONMIN 
optimization  code  [Ref.  2].  Additionally,  guidelines  for 
coupling  the  code  to  an  optimizer  of  the  user's  choice  are 
presented . 


Several  examples  demonstrating  the  program  under  a 
variety  of  conditions  are  presented.  Conclusions  and  recom¬ 
mendations  for  future  work  are  given. 


II.  ANALYSIS 


A.  INTRODUCTION 

When  the  finite  element  method  of  analysis  is  used  to 
design  optimization,  two  objectives  must  be  kept  in  mind. 
First,  the  number  of  analyses  for  the  structure  should  be 
kept  to  a  minimum.  Second,  the  amount  of  gradient  informa¬ 
tion  required  during  the  design  process  should  be  minimized 
to  shorten  run  times  and  minimize  computer  storage 
requirements . 

B.  STATIC  ANALYSIS 

Initial  formulation  of  the  problem  must  include  approxi¬ 
mate  member  areas  in  the  case  of  truss  elements,  and  for 
frame  elements,  characteristic  dimensions  (for  tubular  members 
mean  diameter  and  wall  thickness);  material  properties  (which 
may  be  different  for  each  member) ;  a  set  or  sets  of  external 
loads;  any  non-structural  attached  masses;  and  specified 
joint  support  conditions. 

The  analysis  for  the  stresses  and  deflections  at  the 
joints  must  satisfy  the  conditions  of  equilibrium  of  forces 
at  the  nodes  and  geometric  conditions  of  compatibility  of 
deformation.  In  this  analysis  the  structure  is  assumed  to 
behave  in  a  linearly  elastic  fashion.  The  weight  of  an 
individual  member  is  not  inherently  included  as  part  of  the 


specified  load  conditions;  but,  as  an  option,  half  of  the 
weight  of  each  member  may  be  applied  at  the  member’s  end¬ 
points  as  an  additional  load  in  the  negative  Y-direction. 

For  this  analysis,  the  following  assumptions  are  made: 
truss  and  frame  members  are  treated  as  discrete  entities ; 
truss  elements  have  three  translational  degrees  of  freedom 
at  each  node  and  are  treated  as  pin-connected;  frame  elements 
have  three  translational  and  three  rotational  degrees  of 
freedom  at  each  node  and  are  treated  as  fixed-fixed  beams; 
and  loads  and  reactions  are  applied  at  the  joints  as  shewn 
in  Figure  2.1. 

The  Displacement  (Stiffness)  method  for  finite  element 
analysis  [Ref.  3]  is  utilized  where 

=  £  (3.1) 

and  where  £  is  the  global  stiffness  matrix,  £  is  the  vector 
or  vectors  of  applied  loads,  and  y.  is  the  vector  or  vectors 
of  displacements.  The  method  used  herein  is  an  extension  of 
that  described  by  Felix  and  Vanderplaats  [Ref.  4],  3y 
applying  the  constitutive  stress-displacement  relationships, 
stresses  in  the  elements  may  be  recovered. 

C .  DYNAMIC  ANALYS IS 

When  constraints  on  the  system’s  natural  frequencies 
are  to  be  considered,  the  design  process  requires  the  solu¬ 
tion  of  an  eigenvalue  problem.  The  sub-space  iteration 


Figure  2.1  FORCE  ORIENTATION  CONVENTION  FOR  AN  ELEMENT 


method  of  Bathe  and  Wilson  [Ref.  5]  is  used  to  solve  for  the 
desired  number  of  lowest  eigenvalues  and  the  associated 
eigenvectors.  This  method  is  reasonably  efficient  for  a 
small  number  of  lowest  frequencies  for  a  large  problem,  and 
is  well  suited  for  re-analysis  when  small  changes  are  made 
in  the  design. 

D.  GRADIENTS 

Gradients  are  currently  calculated  with  respect  to  member 
sizes  and  coordinates  only  by  finite  difference  techniques. 
Inclusion  of  the  capability  to  calculate  gradients  analytic¬ 
ally  is  identified  as  a  needed  extension  to  this  work. 


III.  OPTIMIZATION 


*>  »  v  •  *  ,  «  « 


A.  INTRODUCTION 

The  primary  objective  of  structural  optimization  is  to 
design  systems  that  efficiently  perform  specified  purposes. 
Selection  of  a  specific  optimizing  algorithm  must  take  into 
account  the  following:  1)  the  structure  should  be  analyzed 
as  few  times  as  possible;  2)  the  algorithm  should  minimize 
the  amount  of  gradient  information  required;  3)  the  algorithm 
should  provide  reasonable  assurance  that  an  optimum  or  near¬ 
optimum  design  will  be  reached. 


B .  GENERAL  FORMULATION 

The  general  statement  of  inequality  constrained  minimiza¬ 
tion  is  as  follows: 


Minimize 

F(X) 

*«# 

Subject  to: 


G.(p  1  0 


xj  <  x.  <  *J 


j  =  l,m 


i  =  l  ,n 


(4.1) 


(4.2) 


(4.3) 


where  F(X)  is  the  objective  function  to  be  minimized.  The 
functions  Gj(X)  are  the  set  of  inequality  constraints  to 
be  met.  The  vector  X  contains  the  design  variables. 
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The  inequality  constraints,  Gj(X)  £  0  j=l,m,  must  be 
satisfied  for  the  design  to  be  accepted  as  feasible.  Side 
constraints ,  x^  and  x^ ,  are  lower  and  upper  bounds  on  the 
design  variables.  The  objective  function  must  be  minimized 
as  much  as  possible  while  still  satisfying  all  inequality 
constraints.  In  the  case  that  it  is  not  possible  to  satisfy 
all  constraints ,  the  most  nearly  feasible  solution  must  be 
found.  Felix  and  Vanderplaats  [Ref.  4]  is  an  excellent 
source  for  the  basic  structural  design  formulation. 

C.  DESIGN  VARIABLES 

The  vector  X  contains  the  design  variables;  in  this  case, 
member  characteristic  dimensions  for  frame  elements,  cross- 
sectional  areas  for  truss  elements ,  and  spacial  coordinates 
of  the  joints.  The  user  may  elect  to  optimize  the  structure 
weight  with  respect  to  any  combination  of  the  above  variables. 

For  truss  elements,  where  normal  stress  is  dependent  only 
upon  the  magnitude  of  the  cross-sectional  area  of  an  element 
and  not  the  distribution  of  that  area,  it  is  sufficient  to 
use  the  area  as  the  design  variable  as  long  as  the  Euler 
buckling  stress  can  be  related  to  the  cross-sectional  area. 

For  frame  elements  where  stresses  are  dependent  on  area 
distribution,  more  flexibility  is  allowed  by  varying  two 
characteristic  dimensions  and  allowing  the  code  to  calculate 
the  section  properties  from  these. 
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In  the  case  of  tubular  elements,  the  characteristic 
dimensions  are  mean  diameter  and  wall  thickness;  from  which 
area,  bending  and  polar  moments,  and  maximum  radial  dimen¬ 
sions  are  generated. 

The  optimum  geometry  problem  requires  that  the  joint 
coordinates  be  design  variables.  The  x,  y,  and  z  coordinates 
of  a  joint  are  treated  as  separate  design  variables. 

In  many  cases  it  is  desirable  to  link  a  set  of  design 
variable  together  to  preserve  symmetry,  limit  the  number  of 
variables  to  be  solved,  or  to  limit  the  number  of  unique 
elements  to  be  manufactured.  The  code  has  provisions  for 
design  variable  linking  by  which  two  or  more  variables  may 
be  linked  in  equality  or  some  fixed  ratio. 

D.  OBJECTIVE  FUNCTION 

The  objective  function  under  consideration  is  weight 
NE 

F(X)  =  Z  p.A.L.  (4.4) 

i=l  1  1  1 

where  is  the  weight  density  (in  consistent  units)  of  the 
material  of  the  ith  element,  is  the  cross-section  area 
of  the  ith  element,  and  is  the  length  of  the  ith  element 
and  NE  is  the  number  of  elements  in  the  structure. 

E.  CONSTRAINTS 

This  code  is  designed  to  accommodate  constraints  on  maxi¬ 
mum  normal  stress,  Euler  buckling,  local  or  shell  buckling. 


maximum  joint  displacements,  and  minimum  natural  frequencie 
of  the  structure.  All  constraint  values  are  normalized. 


1.  Stress : 

For  truss  elements,  maximum  normal  stress  is 

calculated  directly  from  the  element  tensile  forces.  For 
frame  elements*,  tensile  stress  oA  ,  maximum  bending  and 
shear  stresses  and  are  calculated.  Maximum  and 
minimum  normal  stresses  are 

°max  =  max  magni"tuc^e  of  0A  +  ^cr”A  ^  +  4a^)  (4.5) 

2 

j - 5 - t~ 

amin  =  m^n  magnitu<:^e  °f  1  +  4<7T  ^  (4.6) 

2 

The  upper  and  lower  bounds  on  stress  may  be  different  for 
each  member,  but  are  taken  to  be  the  same  for  every  loading 
condition . 

2 .  Local  or  Thin  Shell  Buckling 

The  stress  at  which  local  or  thin  shell  buckling 
occurs  is  given  by: 


0.4Ei 


(4.7) 


mi  1  'i 


where  the  subscript  i  corresponds  to  the  member  number, 
is  Young's  modulus,  Dm>  is  the  element's  mean  diameter,  and 
t-  is  the  element's  wall  thickness. 


F.  DESIGN  VARIABLE  BOUNDS 

Side  constraints  are  imposed  on  the  design  variables 
as : 

CD*  £  CDi  <  CDj  (4.8) 

•  £  n 

CD  is  the  characteristic  dimension  and  CD^  and  CD^  are 
the  minimum  and  maximum  allowable  characteristic  dimensions 
of  the  ith  element,  and  are  taken  to  be  the  same  for  all 
load  cases.  If,  as  is  the  case  in  a  tubular  element,  geo¬ 
metry  dictates  that  some  relationship  between  the  design 
variables  cannot  be  exceeded,  i.e.,  thickness  cannot  exceed 
the  mean  diameter,  the  user  must  arrange  the  bounds  to  pre¬ 
clude  such  an  occurrence . 

G .  OPTIMUM  GEOMETRY 

Joint  coordinates  are  treated  as  design  variables  with 
a  separate  design  variable  for  each  coordinate  direction 
of  a  node.  Coordinate  design  variables  may  also  be  linked. 
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IV .  PROGRAM  FEATURES 


A.  INTRODUCTION 

Computer  codes  each  have  their  own  features  and  formats 
with  which  the  user  must  become  familiar  if  he  is  to  use 
the  code  easily  and  efficiently.  The  SADX  code  developed 
in  the  course  of  this  research  has  been  designed  with 
various  options  which  are  explained  in  this  chapter. 

Chapter  V  contains  a  User's  Manual  with  sample  data  for  a 
typical  problem  that  might  be  solved  with  this  code:  a 
truss-braced  cantilever  beam.  This  problem  along  with 
other  numerical  examples  will  be  presented  in  detail  with 
results  in  Chapter  VI. 

The  SADX  code  was  written  to  be  used  as  a  stand  alone 
analysis  program  or  as  an  analysis  code  that  could  be  easily 
be  coupled  to  an  optimizer  (of  the  user's  choice)  through 
simple  modifications  to  the  main  driver  program. 

B .  FEATURES 

When  the  user  supplies  member  areas,  section  types, 
characteristic  dimensions,  material  data,  connectivities 
and  joint  coordinates,  along  with  a  set  of  program  control 
parameters;  the  analysis  mode  will  calculate  the  weight  of 
the  structure.  The  addition  of  one  or  more  se^s  of  loading 
parameters  will  result  in  the  calculation  of  resultant  joint 


23 


displacements,  member  stresses,  and/or  forces  for  each  load 
case  along  with  the  desired  number  of  structure  natural 
frequencies  and  modes.  Design  variables  may  be  chosen  as 
truss  element  areas,  frame  element  characteristic  dimen¬ 
sions,  and  joint  coordinates.  When  coupled  to  an  optimizer, 
the  code  will  minimize  the  weight  of  the  structure  and  print 
the  final  optimization  information.  The  code  is  designed  to 
be  as  simple  to  use  as  possible  while  retaining  the  flexi¬ 
bility  for  use  on  a  variety  of  problems . 

The  code's  modular  construction  allows  the  user  to  use 
frame  element  cross-sections  other  than  tubular  elements. 
This  is  done  by  reading  in  two  characteristic  dimensions  for 
each  frame  element  along  with  a  section  type  identifier. 
Subroutine  SADX85  is  called  to  calculate  the  section  proper¬ 
ties,  area,  maximum  radial  dimensions,  and  bending  and  polar 
moments  of  inertia.  The  user  may  augment  this  subroutine  to 
calculate  section  properties  for  whatever  section  type  he 
may  choose  to  work  with. 

Many  computer  codes  require  the  definition  of  an  auxil¬ 
iary  node  to  orient  the  principle  axis  of  a  non-axially 
symmetric  element.  In  this  code  it  is  assumed  that  the 
element's  local  z-axis  is  parallel  to  the  global  x-z  plane. 

Pseudo-dynamic  storage  is  used  to  allow  storage  of  most 
data  in  one  integer  and  one  real  array  for  more  efficient 
use  of  storage. 


Various  print  control  options  are  available  to  tailor 
the  printed  output  to  match  the  user's  desires. 

Design  variable  linking  is  available  to  allow  elements 
and  joints  to  be  grouped  to  maintain  symmetry,  limit  the 
number  of  independent  design  variables,  or  reduce  the 
variety  of  member  sizes  generated. 

Optimization  may  be  performed  with  respect  to  member 
size,  with  respect  to  structure  geometry,  or  with  respect  to 
both. 

Structures  may  be  comprised  of  truss  elements,  frame 
elements,  or  a  mix  of  the  two  types. 

Structures  may  be  optimized  for  multiple  load  cases  with 
constraints  imposed  upon  any  combination  of  maximum  normal 
stress,  maximum  joint  displacements  and  rotations,  Euler 
and  local  buckling,  and  a  specified  number  of  minimum  natural 
frequencies  in  free  vibration.  Separate  displacement  con¬ 
straints  may  be  imposed  for  each  load  case.  Both  consis¬ 
tent  and  lumped  mass  options  are  coded.  Either  forces  or 
stresses  or  both  can  be  output.  The  user  may  decide  whether 
to  include  the  structure's  weight  and  the  fixed  masses  as 
loads  applied  to  the  structure. 

C .  EXAMPLE 

The  following  example  of  the  truss-braced  cantilever 
beam  presented  in  Tables  (I-V)  demonstrates  some  of  the 
options  available  in  the  code. 
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TABLE  I 


TRUSS -BRACED  CANTILEVER  BEAM:  INPUT  CONTROL  PARAMETERS 


INPOT  PARAMETERS  POR  STRUCTURAL  ANALYSIS  AND  DESIGN 
ROUTINE,  "SADX" 

TRUSS-BRACED  CANTILEVER  BEAM:  EXAMPLE  1 
CONTROL  PARAMETERS 

TOTAL  NUMBER  OP  ELEMENTS,  NE=  4 

TOTAL  NUMBER  OP  BAR  ELEMENTS,  NEB*  2 

TOTAL  NUMBER  OP  FRAME  ELEMENTS,  NEF*  2 

TOTAL  NUMBER  OP  JOINTS,  NJ=  5 

JOINT  MAXIMUM  DEGREES  OF  FREEDOM,  JCM*  6 

JOINT  CONSTRAINT  VARIABLE,  NCJ*  3 

NUMBER  OF  MATERIAL  TYPES,  N MT  =  2 

HUMBER  OP  LOAD  CONDITIONS,  NLC*  2 

NUMBER  OF  EIGENVALUES,  NEIG=  1 

NO.  OF  EIGENVALUES  CALCULATED,  NEIG1*  2 

NUMBER  OF  FIXED  MASSES,  NFM ASS*  1 

EULER  BUCKLING  CONSTRAINT  ID,  NEU3C=  1 

LOCAL  BUCKLING  CONSTRAINT  ID,  L8UCK*  1 

NO.  CP  DISPL.  CONSTRAINTS,  ND5PLC*  2 

NO.  OF  PREQ.  CONSTRAINTS,  NFREQ*  1 

LUMPED  MASS  OPTION,  LM ASS*  0 

FORCE/STPESS  PRINT  OPTION,  NSTRES*  2 

OPTIMUM  SIZE/BOTH/GEOM  OPTION  IDVCLC*  2 

STRUCTURE  WEIGHT  AS  LOADS  OPTION  NSTW=  1 

FIXED  MASSES  AS  LOADS  OPTION  NFMW*  1 

ACCELERATION  DUE  TO  GRAVITY,  GRAV  =  0.3864QE+03 

EIGENVALUE  CONVERGENCE  TOLERANCE,  EPSEIG  =  0.100003-03 
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TABLE  II 


TABLE  III 

EXAMPLE  1:  SAMPLE  DISPLACEMENT  AND  FORCE/STRESS  OUTPUT 


»«.  C»  IU  HIPMU  • 


T 

T 

UWO  CCM 

1 

lOAt^CCfcf 

tin  me  fwci  m  i 

-e.?JCJF*04 

0«i94«E*04 

’■“lilsOISL"  * 

•mi  im$s 
•4*H2IT«0) 
0.444M»tft 

x^m 

2 

•4* 

S4lii*&$ 

NC.  W  < 

WUM  CUMillt  • 

2 

lUMKf 

kQAf  COM 

DC Of 

TINS  III  TC»C|  M)«CC 

IN  X  5|XfX  IN  r  0|R€X 

3NEM  FORCI 

In  x  dim* 

TOC  St  CM  N9WNT 
AAOUT 1  AIM 

2SJSI-?  7S?fT  IiffiV‘?!S3NT  ,»2„ 

1 

1 

tit* 

tBW* 

4:t2ta:8i 

-2:ll5ii:st 

8:8 

-8:??:u:2J  8:m?i:82 

*0s1«2Zf*0i 

0.1*C!I»09 

8 

i 

hf&s 

IftKKtt 

1;ti8«48 

2:8 

8:414:3  :2:lU?f4l 

0.U4Tt»09 

o0.2T«0|*09 

1  UNCUT 

iOAC  CCNO 

sect 

IS-l'ifsIS***  BTiV 

INCAN  COCCI 
IN  z  OIMX 

TMMYnr  1SS811  ISt*'  IBSmSB"  sM» 

4 

4 

\ 

hISn 

1:423:22 

14188:8 

*8::?3i*42 

2:2 

m'8  *8:224141 

-SrillSrS? 

uJm'j.s  l»c  t..«S^Te.sg:i?*tf:a 

4:il«!4l 

4:128214: 

2:8 

=2:12?S8:»  .1-2888141 

4.-2*234? 

!5lt»ili8l,<  e.liaaH'Of  cps 

»  -g.jorm-ol  -C.tc7»?I-<il 
s  -4.Jj*3)*-aI  -o.ihsit-ot 


U»*c  CCADITIO*  | 

5»'lRBS?s  «•  '{m 

i  m**  :iiPKi  flw* 


ICAO  CCAOIT1CA 


*CT  AST  X 

0.0 

cJo 

B 


ROT  AST  V 

0*2066-05 

0.206E-05 

0.0 

0.0 


5jr«8!Ks  jSSSff '  w, 

|  '-hllfzSi  rfc'Hltt!  til 


ROT  AST  X 


ROT  AST  T 

0.0 

0.4166—04 
0. 6686-04 
0.0 
Q.O 


ROT  AST  l 
O.Q 

-0 .5676-04 
-0. ST5E-04 
0.0 
0.0 


ROT  AST  I 


04  -0.8476—04 

04  -0.1436-03 

0.0 

0.0 


T  AST  l 

:S 


tmt  | 

*646846-04 

.£4264i-0J 


TABLE  IV 


EXAMPLE  1:  SAMPLE  CONSTRAINT  OUTPUT  INFORMATION 


DISPLACEMENT  CONSTRAINTS  START  AT  NUMBER  2  ANO  STOP  AT  NUMBER  S 

TRUSS  ELEMENT  STRESS  CCNSTRAI NTS 

START  AT  NUMBER  6  ANO  STOP  AT  NUMBER  IT 


FRAME  ELEMENT  STRESS  CONSTRAINTS 

START  AT  NUMBER  18  ANO  STOP  AT  NUMBER  41 


THE  NUMBER  CF  CONSTRAINTS  NCTOT 


41 


i! 

9) 

1 


-.4CE882E+C2 
-.1741116+C1 


H 


__  41111 
-.1012366+C1 
-• 5757226+CC 
-.9C1716E+CC 
211  “•471251 5+CO 
251  -.14C20ll+Cl 
291  —.9556416+00 
23)  -.5559996+00 
37)  -.1959946+01 
41)  -.951795E+C0 


-• 1610556+01 
-.581588E+00 


-.5E764CE+00 
102139E+01 
124755E+C1 
■• 118559E+01 
■.  225C14E+CO 
■•  171444E+01 
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—•581741 i 
-.1327661 


TABLE  V 


EXAMPLE  1:  FINAL  OPTIMIZATION  INFORMATION 


FINAL  OPT  If*  12AT1CN  INFORMATION 


oej 


C.5E75C7E+C3 


DECISION  VARIABLES  (X-V ECTOR) 

C.348Cf£+01  0. 849S2E+01  0.483896 *01 
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•0.123416+01 
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C.45432E+02 

CONSTRAINT  VALUES  (G-VECTGR) 
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71-C. 101846+01  -0.943366+00  ■ 

13 )-C. 975  72e+00  -0.  10214E+01  • 
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371-C. 199956+01  -O.IOOOOE+Ol  « 

THERE  ARE  1  ACTIVE  CONSTRAINTS 

CONSTRAINT  NUMBERS  ARE 

THERE  ARE  0  V1CLATE0  CONSTRAINTS 

THERE  ARE  0  ACTIVE  SIDE  CONSTRAINTS 

TERMINATION  CRITERION 

AeSIOBJI  II-CBJU-1)I  LESS  THAN  OABFUN  FCR  3  ITERATIONS 
NLREER  OF  ITERATIONS  •  12 

OBJECTIVE  FUNCTION  WAS  EVAUIATEO  114  TIMES 

CONSTRAINT  FUNCTICNS  MERE  EVALUATED  114  TINES 

THIS  RUN  RECUIREC  116  STRUCTURAL  ANALYSES 

NUMBER  OF  S6CCNCS  RE GUI  RE 0  FOR  EXECUTION  •  2.79 

WEIGHT  OF  STRUCTURE  GIVEN  AREAS  G  LENGTHS 
WEIGHT-  0.56751E+C3 

TCTAL  WEIGHT  INCLUCING  FIXED  MASSES 
TOTAL  WEIGHT-  C. 83751E+03 


-0. 17411E+01 
-0.10377E+01 
-0.90 172E+00 
— 0. 97270E+CG 
-0.95964E+00 
-0.065566+00 
-0.951806+CO 


JOINT  COCROINATES 
JCINT  X 

1  0.0 

2  0. 1COCOE+C3 

3  O.2C0COE+C3 

4  0.0 

5  O.C 


0.0 

0.0 

0.0 

0. 84259E+02 
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S.o 

.0 


.0 
0.0 
0  .0 
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ELEMENT  INFORMATION  FOR  BAR  ELEMENTS 

ELEMENT- JCINT  RELATIONSHIPS 

ELEMENT  NODE  1  NOCE  2  AREA 

124  C. 3481E+C1 
225  0.84996+01 

ELEMENT  INFORMATION  FOR  FRAME  ELEMENTS 

ELEMENT-JCINT  RELATIONSHIPS 


LENGTH 
0. 1308E+03 
0. 1058E+03 


INC  NCCE1  NCCE2  AREA 


LENGTH 
0.1 0006 +03 
0. 10006+03 


CHAR.OIM. 1 
0 .48396+01 
0.48396+01 


CHAR.0IN.2 

0.21906+00 

0.15586+00 


V.  USER  GUIDE 


A.  INTRODUCTION 

In  developing  any  computer  code  for  engineering  analysis, 
it  is  necessary  to  additionally  develop  concise,  easily 
understood  documentation.  This  SADX  USER’S  GUIDE  is  written 
to  be  easily  understood  by  the  user  having  only  minimal  know¬ 
ledge  of  the  FORTRAN  language .  The  format  follows  that  of 
the  optimization  code,  COPES /CONMIN  [Ref.  2], 

This  chapter  is  devoted  to  acquainting  the  user  with  the 
code  and  necessary  input  data. 

B.  DESIGN  EXAMPLE 

The  simple  example  of  a  four-element  combined  truss /frame 
structure  is  used  to  demonstrate  some  of  the  features  of  the 
SADX  program. 

The  structure  is  shown  in  Figure  5.1,  and  consists  of 
two  tubular  frame  elements  along  the  x-axis  with  two  truss 
braces  to  the  y  and  z  axes  from  the  joint  between  the  frame 
elements.  A  non-structural  fixed  mass  is  attached  at  the 
outboard  end  of  the  second  frame  element  where  two  loads  are 
applied . 

C .  SADX  DATA 

The  SADX  program  reads  data  from  unit  5  and  writes  out¬ 
put  on  unit  6.  Units  30  and  40  are  used  as  scratch  files. 
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OPTIMUM 


4 


The  scratch  file  numbers  may  be  changed  by  changing  two  cards 
at  the  beginning  of  subroutine  SADX01.  The  SADX  program  has 
the  capability  to  read  unformatted  data.  The  following  sec¬ 
tions  define  the  data  which  is  required  by  SADX.  The  data 
is  segmented  into  "BLOCKS"  for  convenience. 

SADX  data  begins  with  a  TITLE  card  and  ends  with  a  END 
card.  Comment  cards  may  be  inserted  anywhere  in  the  SADX 
data  stack  prior  to  the  END  card,  and  are  identified  by  a 
dollar  sign  ($)  in  column  1.  Data  may  be  of  either  the 
"110"  or  "F10.0"  type  or  may  be  free-format  separated  by 
commas  with  no  imbedded  blanks .  Lines  of  formatted  and 
unformatted  data  may  be  intermixed. 

1 .  Formatted  Data 

Formats  are  of  "110"  and  "F10.0"  type.  "I"  formats 
must  be  right  justified,  and  "F"  formats  must  have  the  deci¬ 
mal  point.  The  number  of  cards  read  per  data  block  is  listed 
at  the  bottom  of  each  block. 

2 .  Unformatted  Data 

While  the  USER'S  MANUAL  data  sheet  defines  SADX  data 
in  formatted  fields  of  ten,  the  data  may  actually  be  read  in 
a  simplified  fashion  by  separating  data  by  commas  or  one  or 
more  blanks.  If  more  than  one  number  is  contained  on  an  un¬ 
formatted  data  card,  a  comma  must  appear  somewhere  on  the  card. 
If  exponential  numbers  such  as  2.5+10  are  read  on  an  unfor¬ 
matted  card,  there  must  be  no  imbedded  blanks.  Unformatted 
cards  may  be  intermingled  with  formatted  cards .  Real  numbers 
on  an  unformatted  card  must  have  a  decimal  point. 
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EXAMPLES 


Unformatted  data 
5, 7, 3. 2, 1.3*6, -5,2 
Eauivalent  formatted  data 


5 

7  3.2  1.3+6 

Unformatted 

data 

2 

2,3 

2  3 

Equivalent 

formatted  data 

col  10 

20  30  40 

NOTE:  The  third  line  of  data  contains  no  commas  and  is 
therefore  assumed  to  be  already  formatted. 

Placement  of  more  than  eight  unformatted  data  on  a 
card  will  create  two  (or  more)  formatted  cards  as  required 
Fields  of  zeros  will  be  created  if  more  data  are  required 
than  are  filled  on  an  unformatted  card. 


D.  CONSTRAINTS 

Constraints  are  calculated  and  stored  in  the  G  vector 
as  listed  in  the  following  chart.  The  total  number  of 
constraints 

NCON=NFREQ+2*NDSPLC+NLC*(2*NEB  +  4*NEF  +  NEB  +  2*NEF) 
freq  displ  stress  stress  buckl  buckl 


Where  NCON  is  the  total  number  of  constraints,  NFREQ  is  the 
number  of  frequency  constraints,  NDSPLC  is  the  number  of 
displacement  constraints,  NLC  is  the  number  of  load  cases 
imposed,  NE3  is  the  number  of  bar  or  truss  elements,  and  NEF 
is  the  number  of  frame  elements.  When  any  of  the  constraints 
are  missing  from  the  G  vector,  all  constraints  are  moved  up 
in  the  vector.  For  example,  if  there  is  no  frequency  con¬ 
straint,  then  a  displacement  constraint  would  fill  the  first 
location  of  the  G  vector. 

CONSTRAINTS  ARE  STORED  IN  THE  G  VECTOR  IN 
THE  FOLLOWING  ORDER: 

NFREQ  FREQUENCY  CONSTRAINTS 

2 *NDSPLC  JOINT  DISPLACEMENT  CONSTRAINTS 

STRESS  CONSTRAINTS  ARE  STORED  ELEMENT  BY  ELEMENT 

FOR  A  GIVEN  ELEMENT  CONSTRAINTS  ARE  STORED  BY  LOAD  CASE 

FOR  A  GIVEN  TRUSS  ELEMENT  AND  LOAD  CASE, 

CONSTRAINTS  ARE  STORED: 

TENSILE  STRESS  LOWER  LIMIT 
TENSILE  STRESS  UPPER  LIMIT 
EULER  3UCKLING  STRESS  LIMIT  (IF  APPLICABLE) 

FOR  A  GIVEN  FRAME  ELEMENT  AND  LOAD  CASE, 

CONSTRAINTS  ARE  STORED: 

NORMAL  STRESS  AT  LOW  NODE  LOWER  LIMIT 
NORMAL  STRESS  AT  LOW  NODE  UPPER  LIMIT 
NORMAL  STRESS  AT  HIGH  NODE  LOWER  LIMIT 


NORMAL  STRESS  AT  HIGH  NODE  UPPER  LIMIT 

EULER  BUCKLING  STRESS  LIMIT  (IF  APPLICABLE) 

LOCAL  BUCKLING  STRESS  LIMIT  (IF  APPLICABLE) 

E .  EXAMPLE 

The  initial  configuration  of  the  braced  cantilever  beam 
is  shown  in  Figure  5.1  Stress  constraints  are  imposed  as 
well  as  constraints  on  Euler  and  local  buckling,  displace¬ 
ment,  and  first  fundamental  frequency.  A  non-structural 
fixed  mass  is  applied  at  the  tip  of  the  beam,  and  two  load 
conditions  (PI,  and  P2)  are  imposed.  The  structure's  own 
weight  will  be  considered  as  an  imposed  load  as  will  be  the 
fixed  mass . 

The  linking  of  design  variables  is  demonstrated  by  link¬ 
ing  the  mean  diameters  of  the  two  frame  members  (Dl,  and 
D2  with  D1=D2).  Member  size  variables  are:  truss  member 
areas  (Al,  and  A2) ,  frame  member  mean  diameters  (Dl,  and 
D2),  and  frame  member  thicknesses  (Tl,  and  T2).  Member 
sizing  DESIGN  VARIABLES  are:  XA(1)=A1,  XA(2)=A2,  XA(3)=D1=D2 
XA( 4 ) =T1 ,  XA( 5 ) =T2 . 

Geometry  variables  are  the  attachment  points  (joints  4 
and  5)  of  the  two  truss  members  on  the  y  and  z  axes  (Y4,Z5). 
Coordinate  DESIGN  VARIABLES  are:  XC(1)=Y4  and  XC(2)=Z5.  The 
structure's  weight  is  then  optimized  with  respect  to  member 
sizes  and  structure  geometry  for  a  total  of  seven  design 
variables . 


1 .  Properties /Conditions 

Two  material  types  are  used:  type  1,  aluminum,  is 
used  for  the  truss  members;  and  type  2,  steel,  is  used  for 
the  frame  members.  The  weight  densities  of  the  materials 
(p)  are: 

type  1  =  0.1  lb /in 3 
type  2  =  0.3  lb /in 3 

The  Young's  moduli  of  the  materials  (E)  are: 
type  1  E  =  10.0E+S  psi 
type  2  E  =  29.0E+6  psi 

The  non-structural  fixed  mass  attached  at  the  tip  of  the 
structure  weights  250  lb.  The  applied  loads  are: 

PI  1000.0  lb  in  the  +y  direction 

P2  1000.0  lb  in  the  -z  direction 

The  acceptable  maximum  normal  stresses  are: 
type  1  -25000  psi  <  a  <  25000  psi 

type  2  -36000  psi  <_  <*max  £  36000  psi 
Displacement  limits,  imposed  upon  joint  number  3  (the  tip) 
for  each  load  case  in  the  direction  of  loading,  are: 
load  case  1  y-direction  +/-  3.0  in. 

load  case  2  z-direction  +/-  3.5  in. 

Bounds,  placed  on  the  positions  of  joints  4  and  5  along  the 
y  and  z  axes,  are: 

joint  4  y-coordinate  0.0  inches  to  200.0  inches 

joint  5  z-coordinate  0.0  inches  to  100.0  inches 

The  minimum  natural  frequency  of  the  structure;  is  constrained 
to  be  greater  than  1Hz. 
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2 .  Input  Control  Parameters 


The  following  input  control  parameters  are  given 


ease  of 

following  the 

example : 

NEB  =  2 

NEF=2 

NJ  =  5 

NCJ  =  3 

NMT=2 

IDVCLC=2 

NDJ=2 

NEUBC=1 

LBUCK=1 

NFREQ  =  1 

NFMASS=1 

NEIG=1 

NEIG1=2 

NPRI=0 

NLC  =  2 

NDSPLC=2 

NSTRES  =  2 

NSTW=1 

NFMW=1 

Table  VI  is  a  listing  of  commonly  used  nomenclature. 

The  following  USER'S  MANUAL  is  divided  into  blocks 
A  through  P.  Appearing  directly  below  each  data  field  line 
are  the  parameters  for  the  TRUSS-BRACED  CANTILEVER  BEAM 
example .  It  is  important  to  note  that  the  user  may  choose 
any  consistent  system  of  units . 
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TABLE  VI 


COMMON  VARIABLE  NOMENCLATURE 


A 


CHARDIN  1,2- 

DIR 

E 

EPSEIG 

FC1,...FCN- 

FX 

FY 

FZ 

GRAV 
ID  VC  LC 


IX 

IY 

IZ 

IX  X 

IY  Y 

IZZ 

JN 


L3UCK 


LC 

LH  ASS 


LNO 
LS  5CT 

NCJ 

NDJ 

HD  SGI 
NDSG2 
NDSG3 

ND  S3  4 


MEMBER’S  CROSS -SECTIONAL  AREA 
LOWER  BOUND  ON  DIS PLACEMENTS 
UPPER  BOUND  ON  DISPLACEMENTS 
CHARACTERISTIC  DIMENSIONS  FOR  FRANS  MEMBERS 
FOR  LSECT.EQ.  1:  MEAN  DIAMETER  AND  THICKNESS 
DIRECTION  1=X,  2=1  -  3=Z,  4=rot  about  x, 

5=rot  about  y,  6=rot  about  z. 

YOUNGS  MODULUS 

CONVERGENCE  TOLERANCE  OF  EIGENVALUE 
SOLUTIONS  ( DEFAULT3. 3301) 

LOWER  BOUND  ON  FIRST, SECOND,  ETC. 

SYSTEM  NATURAL  FREQUENCY 
LOAD  FORCES  APPLIED  IN  THE  X  DIRECTION 
LOAD  FORCES  APPLIED  IN  THE  Y  DIRECTION 
LOAD  FORCES  APPLIED  IN  THE  Z  DIRECTION 
ACCELERATION  DUE  TO  GRAVITY  (DEFAULT®  386.4) 
THE  OPTIMIZATION  OPERATION  IDENTIFIER 

1  FOR  OPTIMUM  MEMBER  SIZE  ONLY 

2  FOR  BOTH  OPTIMUM  MEMBER  SIZE  AND  GEOMETRY 

3  FOR  OPTIMUM  GEOMETRY  ONLY 
CONSTRAINT  IDENTIFIER.  I?  NON-ZERO  THE  X-DOF 

IS  CONSTRAINED 

CONSTRAINT  IDENTIFIER.  IF  NON-ZERO  THE  Y-DOF 
IS  CONSTRAINED 

CONSTRAINT  IDENTIFIER.  IF  NON-ZERO  THE  Z-DOF 
IS  CONSTRAINED 

CONSTRAINT  IDENTIFIER.  IF  NON-ZERO  THE 

ROTATION  ABOUT  THE  X-AXIS  IS  CONSTRAINED 
CONSTRAINT  IDENTIFIER.  IF  NON-ZERO  THE 

ROTATION  ABOUT  THE  Y-AXIS  IS  CONSTRAINED 
CONSTRAINT  IDENTIFIER.  IF  NON-ZERO  THE 

ROTATION  ABOUT  THE  Z-AXIS  IS  CONSTRAINED 
JOINT  NUMBER  (GLOBAL) 

KEULER-  EULER  BUCKLING  COEFFICIENT  FOR  BAR 
ELEMENTS 

LOCAL  BUCKLING  CONSTRAINT  IDENTIFIER 

IF  LEUCK.NE.O,  LOCAL  3UCKLIMG  CONSTRAINTS 
WILL  BE  APPLIED  TO  TU3ULAR  FRAME  MEMBERS 
LOAD  CONDITION 
LUMPED  MASS  OPTIONS 

IF  LMASS.EQ.1  CONSISTENT  MASS  MATRIX  USED 
IF  LMASS.NS.O  LUMPED  MASS  MATRIX  USED 
ELEMENT  NUM3S3 

CROSS-SECTION  TYPE  IDENTIFIER 

LSECT.EQ. 1  INDICATES  A  TU3ULAR  MEMBER 
THE  NUMBER  OF  CONSTRAINED  JOINTS 
THE  NUMBER  OF  JOINTS  WITH  (OR  LINKED  TO) 

DESIGN  VARIABLES 

TRUSS  MEMBER  AREA  DESIGN  VARIABLE  NUMBER 
COORDINATE  DESIGN  VARIABLE  NUMBER 
FRAME  MEMBER  CHARACTERISTIC  DIMENSION  1 
DESIGN  VARIABLE  NUMBER 
FRAME  MEMBER  CHARACTERISTIC  DIMENSION  2 
DESIGN  VARIABLE  NUMBER 
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NDSPLC 
NEB 
NEF 
NE IG 
NEIG1 

NEUBC 


NF  NASS 
NFMW 


NFREQ 

NID 

NJ 

NLC 

NLJ 

NMT 

NPRI 


NSTRES 


-  THE  NUMBER  OF  DISPLACEMENT  CONSTRAINT  SETS 

-  THE  NUMBER  OF  TRUSS  ELEMENTS 
*  THE  NUMBER  OF  FRAME  ELEMENTS 

-  NUMBER  OF  PRECISE  EIGENVALUES  TO  BE  FOUND 

-  NUMBER  OF  SIDE  NV  ALU  E/EIGEN  VECTOR  TO  BE 

EVALUATED  DEFAULT  =MIN(2*NEIG,  NEIG+9) 

-  EULER  BUCKLING  CONSTRAINT  IDENTIFIERS. 

IP  NEUBC. NE.O  EULER  BUCKLING  CONSTRAINTS 
WILL  BE  IMPOSED  ON  BAR  AND  FRAME  ELEMENTS 

-  NUMBER  OF  NON- STRUCTURAL  FIXED  MASSES 

ATTACHED  TO  THE  STRUCTURE 

-  FIXED  MASS  WEIGHT  IDENTIFIER 

IFfNFMW.NE.  0)  FIXED  MASSES  WILL  BE 
CONSIDERED  AS  LOADS 


NUMBER  OF  FREQUENCY  CONSTRAINTS 
NUMBER  CF  INDEPENDENT  DEGREES  OF  FREEDOM 
THE  NUMBER  OF  JOINTS 
NUMBER  OF  LOADING  CONDTIONS 

NUMBER  CF  LOADED  JOINTS  FOR  THIS  LOAD  CONDITIO* 
NUMBER  OF  SEPARATE  MATERIAL  TYPES 
INPUT  PRINT  CONTROL 

IF  NPRI. NS. 0  INPUT  VALUES  are  NOT  PRINTED 
IF  NPRI.EQ.5  LOCI/LOCR/IA/RA  ARRAYS  PRINTED 
THE  FOPCE/STRESS  PRINT  IDENTIFIER 

IF  NSTRES.  SQ.O  MEMBER  FORCES  PRINTED 
IP  NSTRES.  EQ.1  MEMBER  STRESSES  PRINTED 
IF  NSTRES.  EQ.  2  BOTH  are  PRINTED 
STRUCTURE  WEIGHT  IDENTIFIER 

IP  (NSTW.  NE.  Of  THE  STRUCTURE'S  OWN  WEIGHT 
WILL  BE  CONSIDERED  AS  LOADS 
POISSON'S  RATIO 
MATERIAL  WEIGHT  DENSITY 
MINIMUM  ALLOWABLE  NORMAL  STRESS 
MAXIMUM  ALLOWABLE  NORMAL  STRESS 
TORSIONAL  MOMENT  APPLIED  ABOUT  THE  X-AXIS 
TORSIONAL  MOMENT  APPLIED  ABOUT  THE  Y-AXIS 
TORSIONAL  MOMENT  APPLIED  ABOUT  THE  2-AXIS 
INITIAL  VALUE  OF  AREA  DESIGN  VARIABLE 
LOWER  BOUNDS  0  N  XA 
UPPER  BOUNDS  D N  XA 

INITIAL  VALUE  OF  COORDINATE  DESIGN  VARIABLE 
LOWER  BOUNDS  D N  XC 


OF  FREEDOM 


CONDITION 


LOADS 


STRESS 

STRESS 

ABOUT  THE  X-AXIS 
A30UT  THE  Y-AXIS 
ABOUT  THE  2-AXIS 


UPPER  BOUNDS  D  N  XC 


DATA  BiOCK  A 


DESCRIPTION:  Title  Card 

Format  and  Example 


TITLE 

FORMAT 

20A4 

TSOSS-BRACED 

CANTILEVER  BEAM 

FIELD 

CONTENTS 

1 


ANY  80  CHARACTER  TITLE  MAY  3E  GIVEN  ON  THIS 
LINE 


DA£A  BLOCK  3 


DESCRIPTION:  Control  Parameters 

I 

|  Format  and  Example 


NEB 

NEF 

NJ 

NCJ 

NET 

IDVCLC 

NDJ 

FORMAT 

7110 

|  2 

2 

5 

3  2  2 

2 

NEUBC 

LBOCK 

NFREQ 

NFMASS 

NSIG 

NEIG1 

NPRI 

format 

7110 

1 

1 

1 

1 

1 

2 

1  ° 

NLC 

NDSPLC 

NSTRES 

NSTW 

NFMW 

FORMAT 

5110 

Ll 

2 

2 

1 

□ 

NOTE:  DEFINITIONS  OF  PROGRAM  INPOT  CONTROL 
PARAMETERS  ARE  LISTED  ON  NEXT  PAGE 


FIELD 


COST  ENT 


B 

i-;. 

i*.1' 

I 

1 


FIRST  C  ARD 

1  NEB-number  of  bar  elements 

2  NEF-number  of  frame  elements 

3  SJ-number  of  joints 

4  NCJ-number  of  constrained  joints 

5  NHT-number  of  seperate  material  types 

6  IDVCLC-design  variable  control  parameter 

If  (IDVCLC.EQ.1)  NDV=RDVAR1 
optimizes  wrt  member  size  only 
If  (IDVCLC.EQ.  2)  NDV=NDVAR1+  RDVAR 2 
optimizes  wrt  member  size  and  geometry 
If  (IDVCLC.EQ. 3)  NDV»ND?AR2 
optimizes  wrt  geometry  only 

7  NDJ -total  linked  and  design  variable  joints 

(i. e.  number  cf  ’movable’  joints) 

SECOND  CARD 

1  NEUBC-Euier  buckling  constraint  identifier 

If  (NEOBC.SE.O)  EULER  buckling  constraints 
will  be  imposed  on  bar  elements 

2  LBUCK-local  buckling  constraint  identifier 

If  (LBUCK. NE. 0)  local  buckling  constraints 
will  be  imposed  on  tubular  frame  elements 

3  NFREQ-number  of  frequency  constraints 

4  NFMASS-nu mber  of  fixed  masses  attached  to 

structure 

5  NEIG-number  of  precise  eigenvalues  to 

be  evaluated 


6  NEIGI-number  of  eigenvalues  to  be  evaluated 

DE? A0LT=l»in.  Of  (2*NEIG  ,  NEIG+3J 

7  NPRI-prirt  control  identifier 

NPRI.ne.O  input  info  not  printed 
NPRI.eq.5  RA/IA/LOCR/LDCI  arrays 
will  be  printed  for  debugging 

THIRD  CARD 

1  NLC-number  of  load  conditions 

2  NDSPLC- number  of  displacement  constraints 

3  NSTBES-fcrce/stress  print  identifier 

If  (NSTRES.EQ .0)  stresses  are  printed 
If  (NSTRES.EQ . 1)  forces  are  printed 
If  (NSTRES.EQ .2)  both  ire  printed 

4  NSTW-structure  weight  identifier 

If  (NSTN.NE.O)  the  structured  weight  is 
considered  as  loads 

5  NFMW-fixed  mass  weight  identifier 

If  (NFMW.NE.O)  the  fixed  masses  are 
considered  as  loads 
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21IA  liOCK  C 


Dynamic  Analysis  Information 


Format  and  Example 


LMASS 

GBAV 

EPSEIG 

FORMAT 

110,2710.0 

0 

386.4 

0.0 

FIELD  CPU  TESTS 

1  LMASS-lumped  mass  option  (if  LSASS.NE.O) 

the  lumped  mass  matrix  is  used. 

2  GRAV-accleration  due  to  gravity 

(default*386.  4  inches/sec*) 

3  EPSEIG-convaroence  tolarar.ee  on  eigenvalue 

solution,  (da  fault®.  POP  1) 


5 


Sill  S12S1  E 


Oeit  this  block  if  NDJ*0  in  BLOCK  B 
DESCRIPTION:  Coordinate  Design  Variable  Linking  Data 

Foreat  and  Exaaple 


JN 

IX 

II 

IZ 

msa 

PCI 

PCZ 

F05NAT 

■ 

4i10,3f 10.0 

m 

Hi 

HI 

1.0 

1.0 

1.0 

□ 

H 

H 

□ 

1.0 

1.0 

1.0 

ZIStD 


CONTENTS 


1 

2 

3 

U 

5 

6 


IX-design  variable  associated 
I Y- design  variable  associated 
IZ-desigr.  variable  associated 
PCX -parti cipati on  coefficient 
?CY-participati cn  coefficient 
PCZ-particioati on  coefficient 


with  x  coord, 
with  y  coord, 
with  z  coord, 
of  x-coord. 
of  y-coord. 
of  z-coord. 


NOTE:  Number  of  cards  raad=MDJ 


J212A  SLock  f 

2IS£SI£TI0N:  Material  Properties 
Fora at  and  Example 


E 

RHO 

SIGMIN 

SIGMAX 

KEULER 

POISSN 

FORMAT 

6F10.0 

■■■■■■ 

■■■■Mi 

■■■■■■■ 

■■■■■■ 

1 1 .  OE+7 

0.  1 

-25000. 

25000. 

2.3 

0.27 

2.9E+7 

Kil 

-360  00. 

3  60  00 . 

2.0 

0.27 

FIELD  CONTENTS 

1  3-roung‘s  Modulus 

2  RHO-material  dansity 

3  SIG  MIN- minimum  allowable  stress 

4  SIGMAX- maximum  allowable  stress 

5  KEtJLSR-Euler  buckling  coefficient 

6  POISSN-Poisson'  s  ratio 


NOTE:  Number  of  cards  read=NMT 


DATA  BLOCK  G 

Omit  this  block  if  NEB=0  in  BLOCK  B 
DESCBIPTIOH :  Bar  Element  Information 
Format  and  Example 


LNC 

NODE1 

NOOE2 

MATCOD 

NDSG1 

A 

LS  ECT 

3 

2 

4 

1 

1 

3.0 

1 

2 

2 

5 

1 

2 

3.0 

1 

FORMAT 
5110 
f 10,110 


CONTENTS 

1  LNO-element  number 

2  NODEI-global  number  associated  wirh  the 

element’s  1st  node 

3  NODE2-glo bal  number  associated  with  the 

element's  2nd  node 

4  MATCOD-mater ial  type  of  this  element 

5  NDSGl-design  variable  number  associated  with 
this  element' s  area 

6  A-member  cross-sectional  area 

7  LSECT-cross-section  type  identifier 

LSSCT.EQ.  1  indicates  a  tubular  member 


NOTE:  Number  of  cards  read*NEB 


PITA  3 LOCK  H 

Omit  this  block  if  NSF=0  in  BLOCK  B 
DESCRIPTION :  Prams  Element  Information  Pormat  and 

Example 


LNO 

NODE1 

NODE  2 

MATCOD 

NDSG3 

NDSG4 

LSECT 

FORMAT 

7110 

3  1  1 

2 

2 

2 

4 

1 

CHARDIM 1 

CHASDIM2 

FORMAT 

2F10 

—— 

1.0 

EDZI 

3 

2 

3 

5 

l_r_ 

1.0 

FIELD  CONTENTS 

1  LNO-element  number 

2  NODEI-global  number  associated  with  the 


element's  1st  node 

3  N0DE2-global  number  associated  with  the 

element's  2nd  node 

4  MATCOD-material  type  of  this  element 

5  NDSG3-design  variable  number  associated  with 

the  element's  1st  characteristic  dimension 

6  NDSG4-design  variable  number  associated  with 

the  element's  2nd  characteristic  dimension 

7  LSECT-crc ss-section  type  identifier 

LSECT.EQ. 1  indicates  a  tubular  member 
NOTE:  Number  of  cards  read  =  NZP 
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DATA  3 LOCK  I 

DESCRIPTION :  Joint  Constraint  Information 
Format  and  Example 


IZZ 

FORMAT 

7110 

FIELD 


CONTENTS 

JN- joint  number 

IX-  x-d is  placement  constraint  identifier. 
IY-  y-displacement  constraint  identifier. 
IZ-  z-dis placea ent  constraint  identifier. 

I XX-  x-axis  rotation  constraint  identifier, 
IXY-  y-axis  rotation  constraint  identifier. 
IZZ-  z-axis  rotation  constraint  identifier, 
if  . NE.  0  -  corresponding  DOF  constrained 


NOTE:  Number  of  cards  read*NCJ 
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JN  I  PX  FI 


0.0  liooo.  I  0.0  I  0.0  I  0„0  I  0.0 


0.0  I  0.0  -1000.  0.0  I  0.0  I  0.0 


FIS12 


CONTENT 


1  NLJ-nuraber  of  leaded  joints  for  t his  load 
condition 

1  JN- joint  number 

2  PX 

3  FT-  Forces  in  the  directions 

4  ?Z 

5  TX 

6  TT-  Moments  about  the  X,T,Z  axes 

7  TZ 

NOTE:  Number  of  cards  read  per  set=NLJ 
Number  of  sets  of  cards  read=NLC 
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DATA  3 LOCK  K 


Omit  this  block  if  NFMASS=Q  in  BLOCK  B 
DESCRIPTION:  Fixed  Mass  Information 

Format  and  Example 


FIELD  CONTENDS 

1  JN- joint  number 

2  HASS-poirt  mass  at  joint  (JN)  in  force  unifs 


NOTE:  Number  of  cards  read=NFMASS 
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D  A£A  BLOCK  L 

Omit  this  block  if  ID7CLC=3 

DESCRIPTION:  Design  Variable  Information 

(MEMBER  SIZE  VARIABLES) 


Format  and  Example 


XA  (1) 

XA(2) 

....  XA(NDVARI) 

FORMAT 

8F10.0 

20.0 

20.0 

mm 

mm 

XAL  (1) 

XAL  (2) 

.  .  .  .  XAL  (NDVAR  1) 

FORMAT 

8F10.0 

f  0*  50 

0.50 

4.0 

0.10 

0.10 

XAU  (1 ' 

XAU  (2) 

.  .  .  .XAU  (NDVAR  1) 

FORMAT 

3F10.0 

25.0 

warn 

25.  0 

2.5 

4.0 

FIELD  CONTENTS 


1  XA-initial  value  of  area  design  variables 

2  XAL-lower  bounds  on  area  design  variables 

3  XAU-upper  bounds  on  area  design  variables 

NOTE:  read  one  value  of  XA,  XAL,  XAU  for  each 
independent  member  size  variable  defined  in 
BLOCKS  3  and  H 


Number  of  cards  read  =  as  reguired 


DATA  BLOCK  M 

Omit  this  block  if  ID7CLC*1 

DESCRIPTION;  Design  7ariable  Information 

(COORDINATE  7ARIA3LES) 


Format  and  Example 


XC  (1) 

XC(2) 

....  XC  (ND7A R 2) 

FORMAT 

8F10.0 

150.0 

50.0 

XCL  (1) 

XCL  (2) 

.  .  .  .XCL  (ND7A 3 2) 

FORMAT 

3F10.0 

|  °*° 

0.0 

XCtJ  (1) 

XCU  (2) 

.  .  .  . XC  0  (ND7A  R2) 

FORMAT 

3F10.0 

2  00.0 

100.0 

FIELD 


CONTENTS 


1  XC- initial  value  cf  coord,  design  variaoles 

2  XCL-lower  bounds  on  coord,  design  variables 

3  XCU-upper  bounds  on  coord,  design  variables 

NOTE;  read  one  value  of  XC,XCL,XCU  for  each 
independent  coord,  variable  defined  in  BLOCK  D 


Number  of  cards  read  =as  required 


DATA  gltOCjC  N 

Oait  this  block  if  NDSPLC=0  in  BLOCK  B 
DESCRIPTION:  Joint  Displacement  Constraint  Information 
Format  and  Example 


KOI 

DIR 

LC 

3L 

3D 

FORMAT 

■ 

3I10,2P10.0 

3 

2 

1 

-3.3 

3.0 

3 

3 

2 

-3.5 

3.5 

FIELD  CONTENTS 

1  JN- joint  number 

2  DIR-direction 

1 =x , 2=y , 3=  z  displacement 
4=x, 5=y ,5=z  rotation 

3  LC-load  condition 

4  BL-lower  bound  on  displacement 

5  BU-upper  bound  on  displacement 


NOTE:  Number  of  cards  read=  NDSPLC 
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DATA  BLOCK  0 


Omit  this  block  if  NFREQ=0  in  BLOCK  3 
DESCRIPTION:  Frequency  Constraint  Information 
Format  and  Example 


I  PCI 

FC2 

FC3 

FOBS AT 

8F10.0 

mi 


FIELD  CDtHJJNIS 

1  FC1 -  lower  bound  on  first  natural  frequency 
constraint  in  Hz.  (cycles  per  second) 

N  FCN-  lower  bound  on  NF3E2-ih  natural 
frequency  constraint  in  Hz. 

(cycles  per  second) 


NOTE:  Number  of  cards  read  *  is  required 


Sill  1£2£S  ? 


DESCRIPTION:  End  card 
Fcraat  and  Exaaple 


FORMAT 

3A1 


MOTE : 


This  card  BUST 
SADX  data. 


appear  at  the  end  of  the 
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VI .  NUMERICAL  EXAMPLES 


A.  INTRODUCTION 

Design  of  three-dimensional  truss  and  frame  structures 
are  presented  herein  and  the  corresponding  numerical  results 
are  summarized  to  demonstrate  the  use  of  the  SADX  code . 

The  examples  begin  with  the  truss-braced  cantiliver 
beam. 

B.  EXAMPLE  1:  TRUSS-3RACED  CANTILIVER  BEAM 

The  simple  truss-braced  cantiliver  beam,  as  shown  in 
Figure  6.1,  has  been  previously  used  for  the  SADX  USER'S 
MANUAL  example.  It  consists  of  two  steel  tubular  frame 
members  with  a  common  outer  diameter  and  different  wall 
thicknesses  arranged  as  a  cantilever  beam  along  the  X-axis. 
There  is  a  fixed  250  pound  mass  at  the  tip  of  the  beam. 

Two  aluminum  truss  members  are  attached  from  the  beam  mid¬ 
point  to  points  on  the  Y  and  Z  axes .  This  structure  is 
designed  for  optimum  member  size  and  geometry  under  a  set 
of  two  load  conditions  and  subject  to  constraints  on  maximum 
stress,  maximum  joint  displacement,  Euler  and  local  buckling 
and  minimum  structure  natural  frequencies.  The  weight  of 
the  non-structural  fixed  mass  and  the  structure's  own  weight 
are  imposed  as  loads.  There  are  five  member  size  and  two 
coordinate  design  variables,  and  a  total  of  41  constraints. 


NAL 


PI 


The  number  of  analyses  required  for  this  design  is  116  using 
2.75  seconds  of  CPU  time  and  terminating  on  the  12th  itera¬ 
tion.  Of  the  analyses  conducted,  84  were  required  for  the 
calculation  of  gradients.  The  weight  of  the  structure, 
including  non-structural  fixed  masses,  is  reduced  from  9542 
pounds  to  838  pounds.  Results  are  given  in  Table  VII. 

C.  EXAMPLE  2:  TWO-TIER  3-D  PORTAL  FRAME  WITH  TRUSS  X-BRACES 
The  two-tier  three-dimensional  portal  frame  with  truss 
member  diagonal  braces,  as  shown  in  Figure  6.2,  is  a  symme¬ 
tric  moderately  sized  case  that  can  be  analyzed  easily  by  a 
variety  of  other  codes.  There  are  20  truss  elements,  16 
frame  elements,  and  four  non-structural  fixed  masses.  The 
material  used  is  steel.  This  structure  is  designed  for 
optimum  member  size  and  geometry  under  a  set  of  three  load 
conditions  and  subject  to  constraints  on  maximum  stress, 
maximum  joint  displacement,  Euler  and  local  buckling,  and 
minimum  structure  natural  frequencies.  The  weight  of  the 
non-structural  fixed  masses  and  the  structure's  own  weight 
are  imposed  as  loads.  There  are  10  member  size  and  five 
coordinate  design  variables  and  a  total  of  493  constraints. 
The  number  of  analyses  required  for  this  design  is  385  using 
120  seconds  of  CPU  time  and  terminating  on  the  21st  itera¬ 
tion.  Of  the  analyses  conducted,  305  were  required  for  the 
calculation  of  gradients.  The  weight  of  the  structure, 
including  non-structural  fixed  masses,  is  reduced  from  9302 
pounds  to  1462  pounds.  Results  are  given  in  Table  VIII. 


TABLE  VII 
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P 


EXAMPLE  1:  FINAL  OPTIMIZATION  INFORMATION 

FINAL  OPTIMIZATION  INFORMATION 
OBJ  •  C .567 5C7E+C3 


THERE  ARE 


1  ACTIVE  CONSTRAINTS 


CONSTRAINT  NUMBERS  APE 

THEPE  ARE  0  VIOLATED  CONSTRAINTS 

THERE  ARE  0  ACTIVE  SIDE  CONSTRAINTS 

TERMINATION  CRITERION 

A E S  (OBJ (  I  l-OB J (  1—1)1  LESS  THAN  DAB FUN  FOR  3  ITERATIONS 
NUMBER  OF  ITERATIONS  ■  12 

OBJECTIVE  FUNCTION  WAS  EVALUATED  114 

CONSTRAINT  FUNCTIONS  WERE  EVALUATED  114 

THIS  RUN  FECUIREC  116  STRUCTURAL  ANALYSES 

NUMBER  OF  SECCNCS  REQUIRED  FOR  EXECUTION  «  2.79 

WEIGHT  OF  STRUCTURE  GIVEN  AREAS  £  LENGTHS 
WEIGHT-  0 .5875 1E+G3 

TOTAL  WEIGHT  INCLUDING  FIXED  MASSES 
TOTAL  WEIGHT-  C.B3751E+03 


TIMES 

TIMES 


JOINT  COORDINATES 
JOINT  X 

0.0 

0.1CCCOE+C3 

3  0 . 2  COCO  E+G3 

4  0.0 

5  0.0 


0.0 

0.0 

0.0 

0  .  84  2  59E+O2 

0.0 


0.0 
0.0 
0.0 
0  .0 

0.4543 2 E+02 


ELEMENT 

1 

2 


ELEMENT  INFORMATION  FOR  BAR  ELEMENTS 
ELEMENT-JOINT  RELATIONSHIPS 
1  NOCE  2 


NOOE 

2 

2 


NOCE 

4 

5 


AREA 

0.  348 1E+-01 
0.8499E+01 


LENGTH 
0.  13C8E+-03 
0. 1098 E +03 


ELEMENT  INFORMATION  FOR  FRAME  ELEMENTS 
ELEMENT-JOINT  RELATIONSHIPS 


LNC  KC0E1 

3  1 

4  2 


NCCE2  AREA 

2  0.332SE+01 

3  0 .2366E+01 


LENGTH 
0. 1000S  +03 
0. 10C0E+03 


CHAR.OIM.l 

0.4839E+01 

0.4639E+01 


CHAR. DIM. 2 
0.2190E+00 
0.1558E+00 


EIGENVALUES  AND  EIGENVECTORS 

MODE  NUMBER 
FREQUENCY 


IIGINVAI^ 
JOINT 


EIGENVECTOR 


1 

2 

•a 

4 

5 


C.41 S8EE+C2  CPS 
C  .2631SE+C3  RPS 

^§grIecgf°freedom 
2 

0.0 

C.  0 

-0.23435E— 01  -C.32985E-01 

-0.50777E-01  -C. 10745 E-01 

-0.23425E-01  -C. 17155E-01 


C.O 

0.0 


0.0 

0.0 

-0.  745C76-01 
0. 158  30E-03 
0. 100  CCE+01 
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D.  EXAMPLE  3:  DD-953  FOREMAST 

Example  three  is  a  redesign  of  the  forward  mast  on  the 
DD-963  of  SPRUANCE  class  destroyer,  Inc  DD-963  foremast  has 
been  chosen  as  the  third  test  case  for  the  following  reasons : 
1)  the  structure  is  typical  of  the  masts  found  on  many  com¬ 
batants  in  the  United  States  and  other  navies,  2)  high  top¬ 
side  weight  reduction  is  desirable  from  a  stability  viewpoint 
for  any  ship,  3)  the  structural  members  are  predominantly 
tubular,  4)  the  problem  is  sufficiently  complex  to  make 
conventional  design  methods  cumbersome,  5)  the  member  size 
and  loading  information  is  available. 

The  structure  as  shown  in  Figures  6.3  through  6.6  is 
constructed  of  172  aluminum  frame  members.  The  material 
used  in  the  structure  is  5036  H32  aluminum,  an  alloy  with 
moderate  strength,  good  weldability,  and  good  corrosion 
resistance . 

This  structure  is  designed  for  optimum  member  size  under 
a  single  load  condition  and  subject  to  constraints  on  maxi¬ 
mum  member  stress,  maximum  joint  displacement,  Euler  and 
local  buckling. 

Some  structural  simplifications  are  made.  The  weights 
of  mast-mounted  radars,  antennas,  and  platforms  are  modeled 
by  a  series  of  fixed  masses  which  are  imposed  as  loads  along 
with  the  structure’s  own  weight.  The  forces  due  to  halyards 
and  wire  antennas  are  applied  as  loads.  Inertial  forces  due 
to  ships  motion  are  calculated  for  the  initial  design  point 
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and  consist  of  the  product  of  a  member’s  half-weight  and  a 
coordinate-dependent  load-factor  applied  at  each  end  of  the 
element  [Ref.  6].  The  load-factor  multiplier  is  dependent 
on  the  distances  from  a  point  at  the  intersection  of  the 
design  waterline  and  the  ship  centerline  amidships  (frame 
264  1/2).  Ship’s  motion  loads  are  applied  for  a  roll  to 
port.  In  that  the  purpose  of  this  example  is  only  to 
attempt  to  solve  a  large  and  messy  problem  rather  than  to 
produce  an  actual  design  no  attempt  is  made  to  include  wind, 
shock,  and  blast  loads.  There  are  34  member  size  design 
variables  and  a  total  of  1054  constraints.  The  number  of 
analyses  required  for  this  design  is  555  using  4482  seconds 
of  CPU  time  and  terminating  on  the  17th  iteration.  Of  the 
analyses  conducted,  544  were  required  for  the  calculation 
of  gradients.  The  weight  of  the  structure,  including  non- 
structural  fixed  masses,  is  increased  from  48,199  pounds 
to  56,746  pounds.  The  structure  as  modeled  was  initially 
infeasible  due,  most  likely  to  the  structural  simplifications 
made  in  the  topmast  area.  Results  are  given  in  Table  IX. 


TABLE  IX 


EXAMPLE  3:  FINAL  OPTIMIZATION  INFORMATION 


FINAL  QPTIMI ZATICN  INFORMATION 

CGJ  *  C.3986C4E+05 

THESE  ARE  2  ACTIVE  CONSTRAINTS 

CONSTRAINT  NUMBERS  ARE 
E60  8(6 

THERE  ARE  0  VIOLATED  CONSTRAINTS 

THERE  ARE  C  ACTIVE  SIOE  CONSTRAINTS 

TERMINATION  CRITERION 

A8Sl0BJm-08J(I-L>  LESS  THAN  CABFUN  FOR  3  ITERATIONS 
NUMBER  CF  ITERATIONS  *  IT 

OBJECTIVE  FUNCTION  HAS  EVALUATED  553  TIMES 

CONSTRAINT  FUNCTIONS  MERE  EVALUATED  553  TIMES 


THIS  RUN  REQUIRED  555  STRUCTURAL  ANALYSES 

WEIGHT  CF  STRUCTURE  GIVEN  AREAS  &  LENGTHS 
WEIGHT*  0  .  39860  £+05 

TCTAL  WEIGHT  INCLUDING  FIXED  MASSES 
TOTAL  WEIGHT*  0.56746E+05 

NUMBER  CF  SECONDS  REQUIRED  FOR  EXECUTION  *  4482.46 

NC.  CF  FRAME  ELEMENTS  *  172 

ELEMENT-JOINT  RELATIONSHIPS 


LNO 

NCCE1 

NGCE2 

AREA 

LENGTH 

CHAR.CIM.l 

CHAR. DIM. 2 

I 

4 

12 

C. 5979E+Q2 

0.10615*03 

C.  19  CO  5*  02 

0.1002 E+01 

2 

■» 

11 

0.  59  79E+C2 

0.1061E+Q3 

C« 19  COE+02 

0.1CC2E+01 

I 

9 

C.5979E+C2 

0. 11 12E*03 

C. 19005*02 

0.10025+01 

4 

2 

10 

C.5S79S+C2 

0. 1112E*03 

C. 1900E+02 

0.10C2E+01 

5 

4 

a 

0.62475+01 

0.1342  E*03 

C.6563E+01 

0. 3030E+00 

6 

3 

6 

0.624 76+ Cl 

0.  12425*03 

C. 65635+01 

0.3030E+00 

7 

1 

a 

0.624  76*0 1 

0. 1369E*0  3 

C. 65635+01 

C.3C3CE+00 

8 

2 

6 

0. 62475*01 

0.1369E+03 

0.65635+01 

0.3030E+00 

9 

a 

12 

C.62475+C1 

0.1295E*03 

C. 65625+01 

0.3030E+00 

10 

6 

11 

0 • 62  4  7E  *01 

0.1295E+03 

C. 05635+01 

0. 303CE+00 

11 

£ 

9 

C ,6247£*01 

0.11325*03 

0.65635+01 

0.3030E+00 

12 

6 

1C 

0.6247E+C1 

C. 1122E+03 

C. 65635+01 

C.3C30E+00 

13 

3 

7 

0.6247E+01 

0.6846  E*0 2 

0.65635+01 

0. 3030F+00 

14 

4 

7 

C.6247E+C1 

0.6846E*02 

C.6563E+01 

0.3030E+00 

15 

7 

11 

0 . 624  76*C1 

0.6846  6*02 

C.6563E+C1 

C. 3C3CE+00 

16 

7 

12 

0.6247E+01 

0.6846  E *0  2 

0.656  3  E+01 

0.3030E+00 

17 

1 

e 

0.62475*01 

0. 1196E+03 

0.65635+01 

0. 3030E+00 

18 

2 

c 

0 .62  4  7E  *C1 

0.1196  E+03 

0.65635+01 

0. 303CE+00 

IS 

5 

9 

0.6247E+C1 

0.1097E+03 

0.65635+01 

0.3030E+00 

20 

5 

1C 

C.624  7E*C1 

0. 1097E+03 

C. 65635+01 

0. 3C3CE+00 

21 

5 

22 

0.6247E+01 

0.21 72  E+03 

0.65635+01 

0 . 2Q30E+00 

22 

1C 

11 

C. 6247E+C1 

0. 21725+03 

C. 65635+01 

0.30206+00 

23 

11 

22 

0.62  4  7E*01 

C. 34C0E+02 

C. 65635+01 

C.3C3CE+00 

24 

e 

1C 

0.62475*01 

0.1894  5*03 

0.65635+01 

0.30306+00 

25 

12 

21 

0. 5979E+C2 

C. 10  £1 E+03 

C. 19C0E+O2 

0. 1002E+0 l 

26 

11 

1£ 

0. 5979E+02 

0.1081 E+03 

C. 19005  +  02 

0. 10C2E+01 

27 

9 

16 

C.  597  SE*C2 

0.1111 E+03 

C. 19005+02 

C. 1002E+0 1 

28 

10 

17 

C.5979E*C2 

0. 1111E+02 

C.  1900E+02 

0. 1 CC2E+01 

29 

9 

21 

0.7358E+C1 

0.2380  E+03 

0.77675+01 

0.3C165+00 

3C 

1C 

IE 

C. 735  86*Ci 

0.2380E+03 

C.7767E+01 

0.3016E+00 

31 

11 

13 

0. 735aE*0l 

0.6846E+02 

C. 77675+01 

0.3  C166+00 

32 

12 

13 

0.7358E*01 

0.6846 E+02 

C.7767E+01 

0 .30166+00 
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33 

13 

IE 

0.7258E+C1 

0.68465+02 

0.7767E+01 

0.3016E+00 

34 

12 

*1 

C. 135EE+C1 

0. 66465*02 

C.  77  6  76  +01 

0. 3  C16E+00 

35 

1C 

16 

0.7358E+01 

0 .2094 E *03 

0.7767E+01 

0.30165+00 

36 

16 

21 

C.  5424E+C1 

0.1 866 E+03 

C. 50245+01 

0.3437E+00 

37 

17 

18 

0.5424E+C1 

C. 18666*03 

C. 5024E+0! 

C.3437E+00 

38 

18 

21 

0.5424E+01 

0.34COE+02 

0.5024E+01 

0.3437E+00 

35 

16 

17 

C. 5424E+C1 

C. 1663  E*03 

C. 50246+01 

0.3437E+00 

40 

12 

20 

0. 5424E+G1 

0. 1222E*03 

C.5024E+01 

0.3437E+00 

41 

11 

15 

C. 54245+01 

0.1232  E*03 

0.50  24E+01 

0.3437E+00 

42 

20 

21 

C.5424E+C1 

0.64566*02 

C.5024E+01 

0.3437E+00 

43 

18 

15 

0.54245+01 

0.6456E+02 

0 .50  24E+01 

0.34375+00 

44 

15 

2C 

C. 5424E+01 

0.84COE+02 

0.50245+01 

0.3437E+00 

45 

21 

25 

0. 5975E+C2 

0. 10  81 E*03 

C.  19CCE+02 

0. 1 OC2E+01 

46 

18 

28 

0.5979E+02 

0.1081 £*03 

C • 19  006+02 

0.1002E+01 

47 

16 

26 

0.5575E+C2 

0. 11125*03 

C. 19005+02 

0.10025+01 

48 

17 

27 

G.5975E+C2 

0.11 12E*03 

C . 19  C0E+02 

0. 1CC2E+01 

45 

16 

25 

G.1058E+C2 

0.2112E+03 

0.9952E+01 

0.3511E+00 

50 

17 

26 

0. 1CSSE+C2 

C. 21125*03 

C.5952E+01 

0.3511E+00 

51 

20 

32 

0. 5424E+01 

0.12095+03 

C.5024E+01 

0.3437E+00 

52 

15 

2C 

C.5424E+C1 

C.  1205  6*03 

0.5024E+01 

0.3437E+00 

53 

20 

31 

0. 105  8E  +  C2 

0. 16266*03 

C.5952E+01 

0.3511E+00 

54 

*  15 

21 

0 • 1058E  *02 

0.16266*03 

C.5952E+01 

0.3511E+00 

55 

18 

24 

0.1058E+02 

0.6846E+02 

C.5952E+C1 

0.3511E+00 

56 

21 

24 

C.10«e5+C2 

0.6846  E*02 

0.9952E+01 

0.3511E+00 

57 

24 

26 

0. 1058E+C2 

0. 6846E*Q2 

C.S952E+01 

0.3511E+00 

58 

24 

25 

0.1058E+02 

0.6846 E*02 

0 .99  52E+01 

0.35116+00 

55 

16 

25 

C. IC58E+02 

0.  1384E+03 

0.9952E+01 

0.3511E+00 

60 

17 

25 

0.10585+02 

0. 13846*03 

C. 59  526+01 

C.3511E+00 

61 

26 

25 

0. 1648E+02 

0.1558  E*03 

0.10825+02 

0 .49486+00 

62 

27 

28 

0. 164EE+C2 

C. 1558E+03 

C. 1082 5+02 

0.4848E+00 

63 

29 

22 

0. 16485+02 

0.1240E*03 

C.1082E+02 

0. 4848E+00 

64 

28 

20 

C. 1648E+02 

0. 1240 E +03 

0.10825+02 

0.4848E+00 

65 

21 

22 

0. 164EE+C2 

0.76845*02 

0. 1082E+02 

0. 4848E+00 

66 

30 

21 

0. 1648E+02 

0.7684  E+02 

0 • 10  82E+02 

0. 4848E+00 

67 

28 

25 

0.164E6+C2 

0.8400E+02 

0.10  82  E+02 

0.4848E+00 

66 

15 

<C 

0.3142E-C6 

a. IQCOE+QO 

C.O 

0.0 

65 

15 

60 

0.31425-06 

0.10006*00 

0.0 

0.0 

7C 

20 

22 

C. 1646E+C2 

C. 96C0E+02 

C.1082E+02 

0.4848E+00 

71 

29 

25 

0. 5975E+G2 

0.96115+02 

C. 19CCE+02 

C. 10C2E+0 1 

72 

28 

27 

0.S975E+C2 

0.9611 E+02 

0.19005+02 

C. 10025+01 

73 

26 

24 

C.39755+C2 

0.98796*02 

C. 19006+02 

0.1002E+01 

74 

27 

36 

0. 5979E+Q2 

0.9879  E+02 

C.1900E+02 

C.10C2E+01 

75 

26 

25 

0. 1667E+C2 

0.1792  E+03 

0.10566+02 

0.5026E+00 

76 

27 

27 

C. 1667E+02 

0. 17926*02 

C.1056E+02 

0. 50266+00 

77 

29 

40 

0.  16675+02 

0.1000 E+03 

0 . 1056E+02 

0.502fcc+00 

78 

<e 

26 

C. 1667E+02 

C. 10C0E+03 

C.1056E+02 

0.5026E+00 

75 

26 

52 

0. 1667E+02 

0.20665+03 

C.  10565+02 

0. 5C26E+00 

80 

27 

53 

0.  1667E  +  C2 

C.2066E+03 

C.1056E+02 

0.5026E+00 

81 

28 

"  a 

0. 1667E+CZ 

0. 63  £2  E+02 

C. 1056E+02 

0.5026E+00 

8£ 

£5 

23 

0. 1667E+02 

0.63825+02 

C. 10565+02 

0.5026E+00 

9* 

*3 

27 

0.  166  7E  +  C2 

0.63826*02 

0.  10  565+02 

0.50265+00 

94 

33 

25 

0. 1667E+C2 

0.63  82  E+02 

C.  10565+  02 

0. 50265+00 

e« 

34 

29 

0. 1667E+02 

0.1289  E+03 

0 . 1056E+02 

0.5C26E+00 

66 

26 

27 

C.  1667E+C2 

C. 1289E+03 

C. 1056E+02 

0.5026E+00 

87 

29 

4  0 

0.  1667E+02 

0.3254E+C2 

C.10  56E+02 

0.50266+00 

88 

37 

23 

0.  1667E  +C2 

0.32545+02 

0.10  56  E+02 

0.5026E+00 

85 

39 

47 

C. 55  7SE+C2 

0. 72  C9E+02 

C.1900E+02 

0. 1002E+0 1 

50 

37 

45 

0.5979E+02 

0.72096+02 

0.1900E+02 

0. 1CC25+01 

51 

24 

42 

0.3575c+C2 

0.7411 E+02 

0. 19  00E+02 

0. 1002E+0 1 

92 

26 

1 2 

C. 3575E+C2 

0.74116+02 

C. 1900E+02 

C.10C2E+01 

93 

24 

47 

0. 16575+02 

0.1446 6*03 

0.1035E+02 

0.5215E+00 

54 

36 

C. 1657E+C2 

0. 14466*03 

C.1035E+02 

0.5215S+00 

95 

40 

48 

0. 1667E+02 

0.74286+02 

C.1056E+02 

0. 5C26E+00 

96 

21 

46 

0. 1667E+02 

0.7428  E+02 

0. 1056E+02 

0.5026E+00 

97 

25 

26 

0. 16675+C2 

C.  71  <6E+02 

C. 1056E+02 

0. 5026E+00 

98 

25 

27 

0. 1667E+02 

0.71 76E+02 

C.1056E+02 

0.50266+00 

99 

26 

^  C 

C. 1667E+C2 

0.1 2 17  S+0  3 

0. 1056E+02 

0. 5026E+00 

100 

27 

?  s 

0.  166  75  *02 

0. 1217  E+03 

C.  1056S+02 

0.5C26E+00 

101 

34 

25 

0. 16676+02 

0.6168  E+02 

0.1056E+02 

0. 50266+00 

102 

1  c 

26 

C. 1667E+02 

0. 6166  E+02 

C. 10565+02 

0.5026E+00 

103 

35 

42 

0. 1667E+02 

0.9138E+02 

C.1056E+02 

0.5C26E+00 

1C4 

35 

43 

0. 1667E+02 

0.9138  E+02 

0 . 10  56E+02 

0.5026E+00 

72 


VII .  CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

An  existing  finite  element  code  was  expanded  to  encom¬ 
pass  the  more  general  case  of  frame  members;  i.e.,  six 
degrees  of  freedom  per  joint.  Combined  truss  and  frame 
structures  were  designed  for  minimum  weight  with  multiple 
load  conditions  considered. 

The  displacement  method  was  used  for  static  analysis 
and  the  subspace  iteration  method  was  used  for  eigenvalues. 

Several  examples  were  considered.  In  every  case  the 
code  worked  as  an  analysis  tool,  and  significant  weight 
reductions  were  obtained  with  the  coupled  optimizer  CONMIN. 

The  SADX  design  code  has  been  shown  to  be  a  useful  tool 
for  ship  mast  optimum  design. 

3 .  RECOMMENDATIONS 

The  following  recommendations  may  be  of  value  for  future 
work . 

1.  The  routines  necessary  to  calculate  gradients  analyt¬ 
ically  should  be  added  to  the  code. 

2 .  The  code  should  be  extended  to  include  plate  and 
shear  elements  and  a  library  of  member  cross-sections. 

3.  An  out  of  core  equation  solver  should  be  added. 
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4.  The  ability  to  specify  multipliers  for  applying  iner¬ 
tial  loads  along  the  three  coordinate  axes  should  be  added 

in  a  fashion  similar  to  that  used  for  applying  structure's 
own  weight  as  loads.  Such  an  addition  would  simplify  dynamic 
load  analysis  and  design. 

5 .  The  method  of  gradient  calculation  should  be  depen¬ 
dent  on  specific  gradients  required  [Ref.  7]  and  [Ref.  8]. 

6.  Gradients  of  frequency  constraints  would  benefit  from 
a  more  efficient  algorithm  [Ref.  9]. 

7.  The  need  for  a  large  scale  public  structural  optimiza' 
tion  code  still  exists . 

8.  The  code  should  be  modified  to  allow  optimum  member 
size  design  followed  by  simultaneous  optimum  member  size  and 
optimum  geometry  design. 
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APPENDIX  A 


DATA  FILES 

|  A.  INTRODUCTION 

This  appendix  contains  the  data  files  used  to  create 
the  test  cases  in  Chapter  VI.  Additionally  the  data  file 
|  for  the  USER’S  guide  in  complete  form  is  presented. 

! 

I 

i 

i 

« 

i 


6 


TABLE  X 


DATA  FILE  TRUSS-BRACED  CANTILEVER  BEAM 


*********1**** **** *2 ******** *3 ****** ***4** ******* 5** *******6*********7********* 

SBLCCKA  TITLE  I  CARO  FORMAT  80H 

1RLS3-BRACEO  CANTILEVER  8E4M:  EXAMPLE  1 

$«**** ***!*** *4*»**2 « ** ****»«3  »***•* **«4 **  •******5*********6«********7********* 


$8 LOCK  a  CCNTPC^  PARAMETERS 


NEB 


|  NEUBC 


LBLCK 


NLC 


NCSPLC 


NJ 


NFREQ 


NSTRES 


NCJ 


NPMASS 


NSTM 


CAROS 

NMT 


NEIG 


FORMAT  3IIO 
I  OVCLC 

NEIGX 


KFMW 


J 

l 


NDJ 


NPRI 


S********l*********2*******»*3****** ***4**  ******* 5*********6********* 7********* 

J9LGCK  C  OYNAMIC  ANALYSIS  INFORMATION  1  CARO  FORMAT  U0.2F10.0 
I  LMASS  GRAY  EPScIG 

$  0  386.1  o! 

S********l**  *******2***’******3******  ***4** *******5 *********6 *********7********* 

SBLCCK  C. JOINT  CCO RO  INATES 


S  NJ  CAPOS 

J  JN  (  X-CCCRO  (  Y-COORO 

1  0.  0. 

2  100.  0. 

3  200.  0. 

4  0.  ISO. 

SB LOCK  E  OESIGN  VARIABLE  LINKING  OAT  A 
J  NO J  CAROS 

*  JN  1  i  1V  l 


Z-COCRO 

0. 
0. 
0. 
0. 
50. 


! 

o 

2 


FORMAT  I 10. 3FI0 


PCX 


FORMAT  4I10.3FIO.0 
PC Y  .  PC2 


1.0  1.1) 


1.0 


1.0 


i.A 

1.0 


j ******** l********* 2 ******** *3 *********4** *******; *********6 *********7********* 

SB LOCK  F  MATERIAL  PROPERTIES  NMT  CAROS  FORMAT  6F10.0 

*  E  RFC  SIGMIN  .  SIGMAX  .  K EULER  .  POISSN 


l.OE*  f 

2. 


LNC 


0.1 

0.3 


I 


I 


-2  5000.  25000.  2!  .2^ 

_  -36000.  36000.  2.  .27 

S*  *******  5,***  ******  2***  ******  3  *********4**  ******* 5*********6*********7********* 

SBLCCK  G  EAR  ELEMENT  INFORMATION  NEB  CAROS  FORMAT  5I10,Fig,I10 

NC0E3  .  MATCOO  NSOGl  AREA  .  LSECT 


NOCE 


I 

4 

5 


3. A 
3.0 


I***** ***!*** ******2 ******** *3  ********* 4**  *******5*********6*********7****** *** 


SBLCCK  H  FRAME  ELEMENT  INFORMATION 

N0CE2  NC0E3  MATCOO 


LNO 


CHAPCIM2 

l.A 

2 


2  *N EF  CAROS  FORMAT  7I10/2F10 
NSOGl  NSDG2  ^  LSECT 

3  4  1 


SC  FARO  IH 1 

*  *.A 

4 

5.0 

s* ******* I*** ****** 2 *********3  ******  ***  4** *******; *********6**** 


SB  LOCK 
S  JN 
$ 


t  JCIfJT  CONSTRAINT  OATA 
IX  IY 


0 

0 

1 


IZ 


NCJ  CAROS 
IXX 

1 

0 

? 

1 


FORMAT  7110 
I  YY 


IZZ 


0 

1 


1* ******* l**« ****** 2 *********3 *********  4**  ******* 5** *******6**** *****7********* 


SBLCCK  J  JCIMT  LOACI NG  CATA 


NLJ 


(CfIT  IF  NLC«C  IN  BLOCK  Bt 

I 
1 


1  CARO 


FORMAT  110 


TABLE  XI 


DATA  FILE  TRUSS-BRACED  CANTILEVER  BEAM  continued 


FV  JL 

iaoo.6 


NLJ  CAROS 
FZ  (  TX 

0. 


FORMAT  HO, 6F10 

J  ^  I  ” 

0.  0. 


3  0.  0.0  -1000.0  0.  0.  0. 

4*444444]*444444442*44*4444*3 44**44 44*4*4  **** 4**5 4* *4 4*4 4*6**4444* **7*4***4**4 

BLOCK  K  FIXED  MASS  CATA  NFMASS  CAROS  FORMAT  U0.F10.0 

ICKIT  IF  N  FM  A  S  $  *  0  IN  BLOCK  9) 

JN  (  MASS 

3  250.0 

4444*4  4  4  1**4  4  4*44  4  2  4*4**444  *3******  4**4**  444441(45  44  444444*6 *44**4**474**444  444 

BLCCK  L  DESIGN  VARIABLE  INFORMATION  CARDS  AS  REQ'O  FORMAT  8F10.0 
(AREA  AND  CIMENSICN  VARIABLES  CMIT  IF  NDVAP1*0J 

XA(IJ  .  .  .  (XA(N0VAR1J  ^  ill 

2C.0  20.0  20.0  2.0  2.o' 

XALIIJ 

.50  • si  4.0  .10  .l&  1  * 

XAU(  I  j 

25.0  35.0  25.0  2.1  4.<J  1  ® 

[4 *4444 44X444 4 4 *4* *2 *** 44444 *3 *444*4 *44 4*4 *4*4 4**5**44444**64*44*444* 7*44*4* *44 

BLOCK  F  DESIGN  VARIABLE  INFORMATION  CARDS  AS  REO'D  FORMAT  9F10.0 


(COORDINATE  VARIABLES 

:cTl)  .  .  .  . XC< NDVAR2 1 


OMIT  IF  NDVAR2-0) 

I  I 


1  I 

4*4 6*444*4 44* 7* 

FORMAT  8F10.0 

I  I 


150.0  50.0 

1  o.i  . 

KZl  .00  .i  ..... 

4*444 44*1 44 4. *4*442 444*4444*3 4444 4* 4*444* ***4***54**4****46*4**44*4474*444444* 

BLCCK  N  JOINT  DISPLACEMENT  CONSTRAINTS  NDSPLC  CAROS  FORMAT  3110, 2F10 
(CHIT  IF  NOSPLC-C  IN  BLOCK  Bl 
JN  t  C IR  i  LC  BL  |  80  I 

3  2  1  -3.0  3.0 

3  3  2  -3.  5  3.5 

44444*441*4444*44424*4*4*4443  *44*4*4444444**4  4**5*4*******6*****4**4  74444*4*44 

BLCCK  C  FREQUENCY  CONSTRAINTS  NO.  OF  CARDS  AS  REO'D  FORMAT  SF10.0 
(CMIT  IF  NFREQ-0  IN  BLCCK  8) 

n  •  . 

1 .0 

********  1*4 **4****2*****4*» *3  ***»•* 4**4**  *444***5**4*444*464*44*4*4*7***4*4  444 

ELCCK  F  ENO  CARD 
NC 


TABLE  XII 

DATA  FILE  TWO-TIER  3-D  PORTAL  FRAME 

$******** i*«*******2'<** *****413  *********  4**  *******5*********6*********7******  *** 

IBLOCK  A  TITLE  1  CARO  FORMAT  80H 

2-TIEF  FRAME  STRUCTURE  WITH  TRUSS  X-8RACSS:  EXAMPLE  2 

j ******** 1 ******** *2 *********3 ****** ***4** **** ***5  ******** *6 *********7 ********* 


I e LOCK  E  CONTRCL 
I  NEB  NEF 


S  NEUBC 

* 


L8LCK 


NDSFLC 


PARAMETERS 

1 1  ^ 

(  NFR6Q 


3  CAROS 
NMT 


NFMASS 


format  alio 
.  IDVCLC 


NEI61 


NSTRSS 


$******** l***  ******2********* 3 *********4**  *******5****»****6*********7********* 

1  ELCCK  C  DYNAMIC  ANALYSIS  INFORMATION  i  CARO  FORMAT  I10.2FIO.O 
S  LMASS  GR/V  EPSE1G 

Sill 
0  2E6.4  0. 

$********!** ****** *2 *********3  *********4** *******5 ********* 6*** ****** 7********* 

SB  LOCK  C  JOINT  COORDINATES 


NJ  CAROS  .  _  _ 

S  JN  X-CCORO  Y-COORO  2-COORO 

1  1  tool  oi  100! 

2  100.  0.  -100. 

3  -100.  0.  -100. 

4  -100.  ,0.  100. 

5  100.  SO.  100. 

6  IO0.  50.  -100. 

7  -100.  50.  -100. 

8  -100.  50.  100. 

9  100.  100.  100. 

10  100.  100.  -100. 

11  -100.  100.  -100. 

12  -100.  100.  100. 

j********l*********2*********3*********4**  * ••* ***5**** *****fc****** ***7 ********* 

S ELCCK  E  CESIGN  VARIABLE  LINKING  DATA 

J  NCJ  CAROS  FORMAT  4U0.3F10.0 

S  JN  17  ,  IY  ,  IZ  ,  PCX  P CY  ,  PCZ 

’  1  4  4  i.4  o.4  1.4 

2  1  0  2  1.0  0.0  -1.0 

3  10  2  -1.0  C.O  -1.0 

4102  -1.0  0.0  1.0 

5  3  4  5  1.0  1.0  1.0 

6  3  4  5  1.0  1.0  -1.0 

7345  -1.0  1.0  -1.0 

834?  -1.0  1.0  1.0 

j****** #*i *********2 ******** *3 *********4******  ***5 ******* *»,.****  ***** 7********* 
IBLOCK  F  MATERIAL  PROPERTIES  NMT  CAROS  FORMAT  6F10.0 

S  E  RFC  SIGMIN  .  SIGMAX  K  EULER  POI  SSN 


100. 

100. 

-100. 

-100. 

100. 

100. 

-100. 

-100. 

100. 

100. 

-100. 

-100. 


Y-COORO 

0. 

0. 

0. 

,0. 

50. 


FORMAT  I10.3FIO 


2.SE+1 


-42C00. 


36  COO • 


|**«*****ir’<********5*********3  ******  **«4**  *******5** ****** *6*********  7********* 

IBLOCK  G  BAR  ELEMENT  INFORMATION  NEB  CAROS  FORMAT  5I10.F10.I10 


IBLOCK  G  BAR  ELEMEN 
I  LNC  NCCE2 
S  I 

2, 5,2,1, 1,5.0, 1 
3, 2, 7,1,1, $.0,1 
4, 6, 3,1, 1,5.0,  1 
5, 3. P.1.1. 5. 0,1 
6, 7, 4,1  ,1,5.0, 1 
7,4,5, 1,1, 5.0,1 
E  » 1,8, 1, 1,5. 0,1 
5, 5, 7,1,1, 5. 0,1 
10,6,8, 1,1, 5. 0,1 
11, 5. 1C.  1,2, 5. C.l 
12,6,5,1,2,5,0,1 
13 ,6  *11, 1,2 , 5. C.l 
14,7,10,1,2  ,5. C, l 


INFORMATION  NEB  CAROS  FORMAT  5110, Flo, 110 

NCDE3  (  MATCOC  NS0G1  AREA  LSECT 

6  1  1  5.0  1 


TABLE  XIII 

DATA  FILE  TWO-TIER  3-D  PORTAL  FRAME  continued 


15. 7. 12. 1.2. 5.  C, l 

16. 8. 11. 1.2.5.  C.l 
17,8.9,1,2.5.0,1 

18. 5.12. 1.2. 5.  C,1 

19. 9. 11. 1.2.5.  C.l 
20 ,10.li,l,2 ,5.0,1 

$44*****4i***44**4^2***4*****3******  *** 4** ******* 5**  ****** *6**** *****7********* 

JBLCCK  h  FRAME  ELEMENT  INFORMATION  2  *NEF  CAROS  FORMAT  7U0/2F10 

i  INC  NCCE2  NCDE3  MATCOC  ,  NSDG1  NS0G2  LSECT 

S  !  I  I  I  f  I  I 

21.1.5.1.3.4.1 

5.C  0.50 

<2,2,6, 1,3,4, 1 

5.0  0.50 

23.3.7.1.3.4.1 

5.0  0.50 

24  ,4, 8, 1,3, 4,1 

5.C  0.50 

25.5.9.1.5.6.1 

5  .0  0.  50 

26, 6, 1C. 1.5, 6.1 

5.0  0.50 

27 .7.11.1.5.6.1 

S.C  C.  50 

28 .8.12.1.5. 6.1 

5.0  *0.50 

25.5.6.1.7.8.1 

.  5.0  0.50 

3C ,6 ,7  ,1 ,7,8,1 


5.C  C. 50 

31.7.8.1.7.8.1 

5.0  0.50 

32. 6. 5. 1.7. 8.1 

5.0  •  0.50 


33, 9, 1C, 1.9, 1C, 1 

S.C  C. SO 

34.10, 11, 1,9, 1C, 1 

S.C  0.50 

35, 11. 12.1,9, 1C, 1 

5.0  0.50 

36, 12, 5, 1,9, 1C, 1 

S.C  C. 5C 

I*  *******  1***  *4  44**2  »**  ******3  *****9  4**4**  »**4*4i*5**  *  ******  6********* 7******  *** 

SB  LOCK  I  JOINT  CCNSTPA INT  OATA  NCJ  CAROS  FORMAT  7110 

J  JN  ,  I*  ,  IY  ,  IZ  ,  I XX  ,  IYY  ,  III 

S  I  I  |  I  I  I  I 

1,1. 1,1, 1,1,1 

2. 1.1. 1.1.1. 1 

3, 1,1, 1,1,1,  l 

4. 1.1.1. 1.1.1 


)••**» 4 4*1*** *4***4 2 *********3******  ***44*  ****  ***5  ********  *6  *********  7**  *****  ** 


S8LCCK  J  JOINT  LCACING  CATA 
I  (CUT  IF  NLC-0  IN  BLCCK  01 

*  NLJ 


5,  1000. 
1C.100C. 
4  , 

<, 500. 
10,500. 


1  CARO 


NLJ  CAROS 

I  TX 


FORMAT  I  10 


FORMAT  I10.6F10 
.  TY  f  TZ 


TA3LE  XIV 


DATA  FILE  TWO-TIER  3-D  PORTAL  FRAME  continued 


5.-5CO. 

e,-5co. 

A  « 

K  .1000. 

S.C.C.C.C.1000. 

12  ,-100C. 

11  9C.OfO.C9— 10C0. 

!•*•****»  1** ****** *2 ******** *3 *********4**  *******5 ******** *6 ********* 7********* 

SBLCCK  K  FIX  SC  MASS  CATA  NFMASS  CARDS  FORMAT  11C,F10.0 

J  ( CUT  IF  NFMASS-0  IN  SLOCK  8) 

*  JN  .  MA<< 

I. .J.0  1  1 

lO.SO.C 

II . 50.0 

12  ,50.  C 

{* *******!*** ****** 2 *********3****** ***4*********5*********6*********7********* 

JeLCCK  L  CESIGN  VARIA8LE  INFORMATION  CAROS  AS  REQ’O  FORMAT  8F10.0 


(AREA  ANC  C IMEN  SI C  N  VARIABLES 
XA<  If  .  .  .  . XA  <  NO VAR  II 


CMIT  IF  NDV AR  1  =  0  ) 


3.5  3.5  5.6  0.5  5.6  0.5  5.6  C.5 

*  xalU°  0,5 

*  0.1  0.1  .lb  0.1  .5&  0.1  . 5 A  0.1 

.50  0.1 

\  XAUIi  I  I  I  III 

3C.0  30.0  30.0  .50  30.0  .50  30.0  G.5 

2C.0  . 5C 

J* ******* 7** ****** *j ******** *3 *********4**  ******* 5** ******* 6********* 7********* 

SB  LOCK  M  OESIGN  VARIABLE  I NFQRNATtCN  CAROS  AS  REQ’O  FORMAT  BFiO.O 


(COORDINATE  VARIABLES 


IOC. 6 
XCLdj 

2C.G 

XCUdj 

16C.0 


100.  A 
20.(1 


XC1N0VAR21 


CMIT  IF  NDVAR2-0I 


ico.i 

20.  A 


!*•*** ***X*« ****** *2 *********3 ********* 4** *******5  ** ******»ft ********* 7 ********* 

SBICCK  N  JOINT  CIS  PLACEMENT  CONSTRAINTS  NOSPLC  CARDS  FORMAT  3110, 2F10 
S  (CUT  ,F  NCSPLC-0  IN  BLOCK  3) 


]*****••* 1 ******** *2 ******** *3 ********* 4** ****  ***5 ********* ^ ******* **7 ********* 

seceCK  C  F7  cQUENCY  CCNSTR A INTS  NO.  OF  CARDS  AS  REQ’O  PORMAT  8F10.0 
S  ( CMIT  IF  NFREC  *0  IN  SLCCK  Bi 

i  “I  *  ‘"V  I  I  «  I  I 

l.c  5.6  10.6 

4********l*********i ******** *3****** ***4** *******5**** *****^*** ****** 7********* 

SELCCK  f  ENO  CARO 
ENC 


TA3LE  XV 


DATA  FILE  DD-963  FOREMAST 


J8L0CK  A  TITLE 

CC-963  CLASS  CE5TROYER  FORWARD  MAST  REDESIGN:  EXAMPLE  l 
S  THIS  IS  A  SIMPLIFICATION  OF  THE  ACTUAL  STRUCTURE  AND  LOADS  8ASEO  ON 
J  NAVSHIPS  DRAWING  NUMBER  128-4535510  AND  NAVSHIPS  SKETCH  NUMBER 
1  80064-128-SK4537579 

J  THIS  ANALYSIS  RECUPES  A  THE  X-AXIS  FORWARD,  Y-AXIS  UP,  AND  2-AXIS 
*  TO  STeD;  A  DIFFERENT  ORIENTATION  THAN  IN  THE  ABOVE  DRAWINGS 


JELCCK  E  CCNT s  CL  PARAMETERS 


IDVCLC 


$  NEUBC 

i 


LBLCK 


NDSFLC 


NFREQ 


NSTRES 


NFMASS 


NEIG1 


J*** *****!** ******* 2 ******** *3  ******  ***4«*  **** *** 5**** ****•£**** ***** 7**« 

JELCCK  C  CYNAM  IC  ANALYSIS  INFORMATION 
Q  386.  A  0. 

SB  LOCK  C  JOINT  COORDINATES 


1 

-240.0 

0.0 

-106.0 

2 

-240.0 

0.0 

1C6.0 

3 

O.C 

0.0 

42.0 

4 

0  .0 

0.0 

-42.0 

5 

-223. C 

54.0 

0.0 

6 

-120. C 

54.0 

68.  4 

7 

-2.64 

54.0 

0.0 

8 

-120. C 

54.0 

-68.4 

9 

-216. C 

1C3.0 

-54.  7 

10 

-216 .0 

103.0 

54.7 

11 

-5.  28 

1G3.0 

42.0 

12 

-5.23 

1C3.0 

-42.0 

13 

-7.92 

162.0 

0.0 

14 

-172. C 

444.0 

0.  0 

5  IS 

AN  unconnected  fixed 

CUMMY 

15  0.1  0.0  0.0 

It  - 152 • 6  216.0  -83.16 

17  -152.6  216.0  83.16 

18  -1C. 5t  216.0  42.0 

19  54.0  216.0  42.0 

20  54.0  216.0  -42.0 

II  -1C.56  216.0  -42.0 

22  -172  . C  736.0  -  13.0 

23  -172.0  736.0  18.0 

24  -13.2  270.0  0.0 

25  -168.8  324.0  0.0 

26  -163.6  32*. 0  -71.76 

27  -168.8  324.0  71.76 

20  -15.84  324.0  42.0 

25  -15.  84  324.0  -42.0 

20  1C8.0  324.0  48.0 

31  168.0  324.0  0.0 

32  1C8.C  224.0  -48.0 

23  -18.12  372.0  0.0 

34  -147.8  420.0  -61.68 

2  5  -1*7.8  4  20.0  0.0 

26  -147.8  420.0  61.63 

27  -20.4  420.0  42.0 

39  12. C  420.0  29.0 

25  -20.4  420.0  -42.0 

4C  12.0  420.0  -29.0 

41  -22.2  456.0  0.0 

42  -132.0  492.0  -54.0 

43  -132. C  452.0  54.0 

44  -172  .0  492  .0  -54.0 

45  -24. C  452.0  42.0 

46  30.0  452.0  36.0 

47  -24.0  492.0  -42.0 

46  20. C  492.0  -36.0 

J  FNC  NCCES  FCF  THE  LGh  ER  QUAORAPOD  CF  THE  FORWARD  MAST 
S  START  THE  UPFER  FCFEKAST  NODES 


32 


'ABLE  XVI 


DATA 

FILE  DD- 

963  FOREMAST  continued 

49 

-172. C 

452.0 

-120.0 

50 

-172. C 

452.0 

120.0 

51 

-228.0 

408.0 

0.0 

52 

-268. C 

452.0 

-56.0 

53 

-2E8.C 

492.0 

56.0 

54 

-240.0 

492  .0 

0.0 

cc 

-172. C 

452.0 

-2  76.  0 

56 

-172.0 

452.0 

276.0 

57 

-213.0 

600.0 

-24.0 

58 

-218.0 

6C0.0 

24.  0 

59 

-172.0 

492  .0 

54.0 

6C  IS 

AN  UNCONNECTED  FIXED  CUMMY  NOCE 

60 

0.0 

0.0 

0.0 

61 

-172.0 

600.0 

0.0 

62 

-172.0 

452.0 

0.0 

43 

-172.  C 

664.0 

0.0 

64 

-172.0 

736.0 

0.0 

45 

-172.C 

856.0 

0.  0 

66 

-172  .0 

576.0 

0.0 

67 

-172.0 

1056.0 

0.0 

66 

-172. C 

736.0 

-2  64.  0 

69 

-172.0 

736.0 

2  64.0 

7C 

-172. C 

736.0 

-156.0 

71 

-172.0 

736.0 

156.0 

72 

-172 .0 

736.0 

-228.0 

73 

-172.C 

736.0 

228.0 
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DATA  FILE  DD-963  FOREMAST  continued 


5.i!  .  25 

5  iff  .11  71  1  29  1 

1  SPSS 5  PLATFORM 

1(0  63  22  1  29  30  1 

5.75  .25 

5 161  63  23  1  29  30  1 

.HI  *11  23  1  29  30  1 

5.75  . 2 . 

163  64  23  1  29  30  1 

5.75  .25 

g  164  64  22  l  29  30  1 

i  loi«eb‘>Iroarp  ANO  EPACES 

165  44  62  1  15  16  1 

7*6^5  «3 75 

7  ,3?  62  1  15  16  l 

7#£*5  #375 

167  49  44  1  31  32  1 

11.50  .50 

163  59  SO  1  31  32  1 

11.50  .50 

165  45  55  1  31  32  1 

11.50  .50 

170  5  0  56  1  31  32  1 

11.  SC  .SC 

.171  14  49  1  33  34  1 

5.75  .25 

177  14  50  1  33  34  l 

5.75  .25 

j*»44**44l4**4444«42**44444443******4**4*********5*********b*********7*** 


18 LOCK  I  JOINT  CONSTRAINT  OATA 


NCJ  CAROS 
j  I XX 


FORMAT  7110 
I  YY 


TABLE  XXII 


DATA  FILE  DD-963  FOREMAST  continued 


38 

35 

40 

a 

43 

44 

45 

46 

47 

48 

49 
JO 


•3 

54 

55 

56 

57 

58 
55 
60 

a 

64 

65 

47 

48 
69 
7C 

71 

72 

73 


0 

0 

0 

0 

c 

0 

0 

0 

0 

0 

c 

0 

c 

0 

0 

c 

0 

0 

c 

0 

0 

c 

1 

0 

c 

0 

0 

0 

0 

0 

0 

0 

a 

0 

a 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

o 

0 

0 

Q 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

c 

0 

0 

0 

0 


BICCK  J  JOINT  LOADING  CATA 

(CUT  IF  NLC-C  IN  8LCCK  81 

NLJ 

76 

JN  FX  FY 

►ALYARC.  AND  h  IRE  ANTENNA  LCAlOS 


FZ 


1  CARD 


NLJ  CARDS 
TX 


FORNAT  110 


FORMATI10.6F10 

TY 


0 

0 

0 

0 

0 

0 

0 

0 

c 

0 

n 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

a 

0 

0 

0 

0 


4 *444 4 441»44 4 444*42 4 *4 4444443 4444444444** 444444454444*4*446444444444 7*«* 


TZ 


JN 


45  -5C. 

-100. 0 

O.C 

0. 

0. 

0. 

50  -50. 

-100.0 

0.0 

0. 

0. 

0. 

*5  -50. 

-100.0 

0.0 

0. 

0. 

0. 

56  -50. 

-1C0.0 

0.0 

0. 

0. 

0. 

67  0. 

FOTICN  LOADS 

-50.0 

0.0 

0. 

0. 

c. 

FX 


3 

4 

5 

6 
7 

a 

5 

10 

II 

\t 

18 

19 


oA 

o.c 

0.0 

o.c 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

o.c 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


FY 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


FZ  I 

160.5 
163.  3 

96.4 

137.5 
1C8.  2 
184.3 

66.  7 

184.3 
414.9 
447.0 

381.6 
381.6 

70.2 

325.4 
5  65.  7 
528.  5 
429.1 
2  26.5 


TX 


I 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

c.o 

0.0 

0.0 

0.0 

0.0 

0.0 


I 

0.0 

c.o 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

c.o 

0.0 

0.0 

0.0 

0.0 

0.0 


TZ 


O.J 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

o.o 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


39 


TABLE  XXIII 


DATA  FILE  DD-963  FOREMAST  continued 


20 

21 

22 

ll 

25 

26 
21 
28 

29 

30 


23 

24 

25 

36 

37 
28 
29 

40 

41 

42 

43 

44 

45 

46 

47 
46 

49 

50 

11 

54 

55 

56 

57 

58 

59 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

i 


81CCK  K  PTXEC  MASS  CATA 

(CMIT  IF  NFM AS  S*C  IK 
JN  MASS 

PLATFOPM  86LCW 


0.0 

0.0 

214.8 

0.0 

0.0 

468.9 

0.0 

0.0 

5  39.0 

OtO 

0.0 

5  29.0 

o.c 

0.0 

?1.6 

o.c 

0.0 

4  74.  7 

0.0 

0.0 

3234.2 

0.0 

0.0 

2296.2 

0.0 

0.0 

18  74.0 
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0.0 

0.0 

0.0 
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0.0 

0.0 
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22  3  1  -6.0  6.0 
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APPENDIX  B 


PROGRAM  ORGANIZATION 


A.  DESCRIPTION 

The  program  organization  is  layed  out  in  the  following 
flow  charts.  The  main  driver  program  CSADXM)  calls  a 
subdriver  (SADXSD),  and  the  optimizer  of  the  user's  choice. 
All  changes  required  for  replacement  of  the  optimizer  are 
made  in  SADXM,  thus  allowing  for  easy  testing  of  several 
optimizers  on  the  same  problem. 

SADXSD  may  be  called  from  the  main  for  input,  analysis, 
and  output.  Printed  output  may  vary  as  the  user  requires. 

A  complete  listing  of  all  subroutines  and  their  functions  is 
given  in  Table  XXIV. 


TABLE  XXIV 


SUBROUTINE  DIRECTORY 

SXT5XH -  DRIVER  PSOGRXIT'FOiruSING  THE  ABOVE 

SUBROUTINES.  S&DXM  MAI  BE  COUPLED 
TO  OPTIMIZES  OF  USER'S  CHOICE. 

SADXSD  SUBDRIVER  PROGRAM  FOR  COUPLING 

SADX  ROUTINES  TO  SADXM 

SADX01  THIS  ROUTINE  READS  AND  PRINTS  INPUT  DATA 

AND  ORGANIZES  PSEUDO-DYNAMIC  STORAGE 
ALLOCATION 

SADX02  BUILDS  VECTORS  JC  AND  IIR  FOR  FINITE  ELEMENT 

STRUCTURAL  ANALYSIS 

SADX03  BUILDS  THE  12x12  ELEMENT  STIFFNESS  MATRIX 

SADX05  SUPERIMPOSES  THE  ELEMENT  STIFFNESS  MATRIX  EK 

(OR  ELEMENT  MASS  MATRIX  EM)  ON  THE  COMPACTED 
GLOBAL  STIFFNESS  MATRIX  AK  (OR  THE  GLOBAL 
COMPACTED  MASS  MATRIX  AM) 

SADX06  BUILDS  GLOBAL  LUMPED  MASS  MATRIX 

SADX07  LU  DECOMPOSES  SYMMETRIC, POSITIVE-DEFINIT ” 

SPARCE  MATRICES, THE  UPPER  TRIANGLE  OF  WHlv.u 
IS  STORED  r N  MATRIX  AK  (OR  AM)  WITH  LEADING 
ZEROES  NOT  STORED 

SADX08  FORWARD  AND  BACK  SUBSTITUTES  TO  YIELD  A 

SOLUTION  A  SET  OF  LINEAR  EQUATIONS 
(DECOMPOSED  BY  SADX07  OS  EQUIVALENT) 

SADX09  PRINTS  ALL  JOINT  DISPLACEMENTS  FOR  EACH 

LOAD  CONDITION  OF  A  FINITE  ELEM.  STRUCTURE 

SADX11  ROUTINE  TO  ORGANIZE  ANALYSIS 

SADX  1 5  CALCULATES  VALUES  FOR  ALL  DESIGN  AND 

BEHAVIORIAL  CONSTRAINTS  AS  DEFINED  BY 
"SADX”  PROGRAM 

SADX16  CALCULATES  STRESS  IN  TRUSS  ELEMENT  LNO 

UNDER  LOAD  CONDITION  JJ 

SADX 17  PRINTS  STRESSES  AND/OR  FORCES  FOR 

TRUSS  ELEMENTS 

SADX  19  ADDS  ELEMENT  MASS  MATRIX  AA  OF  ELEMENT 

LNO  TO  GLOBAL  MASS  MATRIX  AM  TO  BUILD 
THE  LUMPED  M  ASS  MATRIX 

SADX23  CALCULATES  WEIGHT  OF  TRUSS/FRAME  STRUCTURE 

OR  CALCULATE  WEIGHT  OF  INDIVIDUAL  MEMBERS 

SADX36  CALCULATES  (  XEIG-T*AM*XEIG)  FOR 

GRADIENT  CALCULATIONS  IN  FREQUENCY 
CONSTRAINTS 

SADX37  CALCULATES  EIGENVALUE  GRADIENT  INFORMATION 

IN  FINITE  ELEMENT  STRUCTURAL  ANALYSIS  AND 
DESIGN 


SADX46  READS  INPUT  INFORMATION  FOR  TRUSS 

ELEMENTS 

SADX4  7  TRANSFORMS  THE  ELEMENT  STIFFNESS  MATRIX 

EK  (OR  ELEMENT  MASS  MATRIX  EM)  FROM  LOCAL 
TO  GLOBAL  COORDINATES 

SADX49  SOLVES  REAL  EIGENVALUE  PROBLEMS  USING 

THE  SUES? ACE  ITERATION  METHOD 


SABX50 

3ADX53 
SADX6  2 

SADX71 

SADX72 

SADX78 

SADX80 

SADX81 

5ADX82 

SADX83 

SADX84 

SADX85 

SADX86 


BUILDS  INITIAL  SET  OF  3ASI5  VECTORS 
FOR  EIGENSOL  UTION  31  SEDUCED  BASIS 
METHOD 

PRINTS  MEM3ER  INFORMATION  FOR  TRUSS 
ELEMENTS 

PRINTS  NSIG  EIGENVALUES  STORED  IS 
EIG7AL,  AND  THEIR  CORSES  PONDING 
EIGENVECTORS  STORED  IN  XSIG 

PRINTS  G  VECTOR  OF  CONSTRAINTS 

READS  IN  FRAME  ELEMENT  INPUT  DATA 

BUILDS  3x3  TRANSFORMATION  ARRAY  TRFORM 
FOR  TRANSFORMING  FROM  LOCAL  TO  GLOBAL 
COORDINATES 


CALLS  SADX03  TO  BUILD  THE  ELEMENT  STIFFNESS 
MATRIX;  CALLS  SADX47  TO  TRANSFORM  THE  MATRIX 
AND  CALLS  SADX05  TO  BUILD  THE  REDUCED  GLOBAL 
STIFFNESS  MATRIX 


CALCULATES  STRESS  IN  FRAME  ELEMENT  LNC  UNDER 
LOAD  CONDITION  JJ  (IF  JJ=0  STRESSES  CALCUL¬ 
ATED  FOR  ALL  LOAD  CONDITIONS) 

READS  INPUT  DATA  FOR  FRAME  ELEMENTS  NTTH 
SPECIFIED  SECTION  TYPES 


PFINTS  STRESSES  AND/OR  FORCES  FOR  FRAME 
ELEMENTS 

PRINTS  MEMBER  INFORMATION  FOR  FRAME  ELEMENTS 


CALCULATES  SECTION  PROPERTIES  FOR  FRAME 
ELEMENTS  OF  A  SECTION  TYPE  GIVEN  BY  LSECT 


CALLS  Eli  HER  SADX06  TO  BUILD  THE  LUMPED  MASS 
MATRIX  OR  3UILDS  THE  CONSISTENT  MASS  MATRIX 
BY  CALLING  SADX87  TO  BUILD  THE  ELEMENT  MASS 
MATRIX,  S  AD  X  78  TO  BUILD  THE  TRANSFORMATION 
MATRIX,  SADX47  TO  TRANSFORM  THE  ELEMENT  MASS 
MATRIX,  AND  SADX88  TO  ASSEMBLE  THE  COMPACTED 
GLOEAL  MASS  MATRIX 


SADX87 


CALLS  BUILDS  THE  ELEMENT  CONSISTENT  MASS 
MATRIX 


SADX88 


SADX89 

SADX90 

SADX91 

SETIME 

GETIME 


COW VESTS  UNFORMATTED  DATA  TO  FORMATTED  DATA 
IN  FIELDS  OF  10,  EACH  RIGHT  JUSTIFIED  AND 
ACCEPTS  SCOMMENT  CARDS  IN  DATA 

SOLVES  EIGENVALUE  PROBLEM  |A-ALAMBD*6|  7=0 

SOLVES  EIGENVALUE  PROBLEM 


SOLVES  EIGENVALUE  PROBLEM 

STARTS  EXECUTION  TIMER 
(NGNIMSL  LIBRARY) 

STOPS  EXECUTION  TIMER 
(NONIMSL  LIBRARY) 


UNCLASSIFIED 


OPTIMIZATION  OF  THREE  DIHENSIONRL  COMBINED  TRUSS/FRAME 
STRUCTURESCU)  NAVAL  POSTGRADUATE  SCHOOL  MONTEREV  CA 
G  L  BENDER  OCT  82 

.  F/G  13/12 


TABLE  XXV 


SAD  XU 


PROGRAM  BLOCK  DIAGRAM 


•SBTIHE* 


’SADXS1 


(NONIHSL  LIBRARY) 
SADX88 

‘SAD  XO  V 


rr-  SADXce 
SADXOI-f—  SADX46 

V  SADX72— 
^SADX02 


•SADX85 


SADX85 


•31DX03 


■S  ADX1 


SADX80^SADX78 
\-SADX47 
>*SADX05 
ySADX06- 
yfc-SADX87 
I/— SADX86^-S  ADX78 
IT  V-SADX47 

emn„Crt  VSADX05 
SADXSO 
SADX07 
SADX23 
SADX08 

r— SADX08 
SADX49^—SADX89i 


•SADX19 


■SADX23 


‘SADX90 

•SADX91 


'SADX15> 


•SADX16- 
■  SADX91i 


•S  ADX53 
■SADX84 
•SADXOS 
•SADX79 
•SADX17- 


•SADX78 
►S  ADX03 
*3  ADX78 
•SADX03 


■SADXI61 


■SADX83— SADX81- 


,L-SADX62 

Lrqmilggl}0171 


•SADX78 
•3ADX03 
•S  ADX78 
•5ADX03 
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